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Abstract 

We derive the differential equation governing the asymptotic quasi- 
stationary states of the two dimensional plasma immersed in a strong 
confining magnetic field and of the planetary atmosphere. These two 
systems are related by the property that there is an intrinsic constant 
length: the Larmor radius and respectively the Rossby radius and a 
condensate of the vorticity field in the unperturbed state related to 
the cyclotronic gyration and respectively to the Coriolis frequency. 
Although the closest physical model is the Charney-Hasegawa-Mima 
(CHM) equation, our model is more general and is related to the 
system consisting of a discrete set of point-like vortices interacting 
in plane by a short range potential. A field-theoretical formalism is 
developed for describing the continuous version of this system. The 
action functional can be written in the Bogomolnyi form (emphasiz- 
ing the role of Self-Duality of the asymptotic states) but the minimum 
energy is no more topological and the asymptotic structures appear 
to be non-stationary, which is a major difference with respect to tra- 
ditional topological vortex solutions. Versions of this field theory are 
discussed and we find arguments in favor of a particular form of the 
equation. We comment upon the significant difference between the 
CHM fiuid/plasma and the Euler fluid and respectively the Abelian- 
Higgs vortex models. 
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1 Introduction 

The instabiUties of a plasma embedded in a confining magnetic field [e.g. ex- 
perimental fusion devices, like tokamak) evolve in a geometry that is strongly 
anisotropic. The motion of the electrons along the magnetic field lines is in 
many cases sufficiently fast to produce a density perturbation that has the 
Boltzmann distribution in the electric potential. In these two dimen- 

sional approximation (the dynamical equations are written in the plane that 
is transversal to the magnetic field) may be satisfactory. The Boltzmann 
electron distribution suppresses the convective nonlinearity, leaving the ion 
polarization drift as the essential nonlinearity. This is of high differential 
degree (therefore it is enhanced at small spatial scales) and basically de- 
scribes the advection of the fluctuating vorticity by the velocity fluctuations. 
The differential equation for the electrostatic potential has been derived by 
Hasegawa and Mima \T\. A similar situation appears in the physics of the at- 
mosphere and the differential equation for the streamfunction of the velocity 
field has been derived by Charney [2]. 

Many plasma instabilities (in particular in tokamak) depend on this non- 
linear term and will likely show some common aspects. In numerical simu- 
lations of the Charney-Hasegawa-Mima (CHM) equation it has been proved 
that the plasma is evolving at large times (in the absence of external driving 
forces and starting from an irregular flow pattern) to states that are char- 
acterised by a regular form of the potential. The evolution toward a very 
regular pattern of vortical flow is also characteristic to the incompressible 
ideal fluid, described by the Euler equation [3], [1], [S], [0], [7], [Bj, |9j. In 
this case the time-asymptotic states consists of few vortices, of regular shape, 
with very slow motion. The streamfunction obeys, in these states, the sinh- 
Poisson equation. In the case of CHM equation, the cuasi- stationary states 
also consist of structures but it has not been possible to derive an equation 
for the streamfunction [TO]. There are seveal studies of the CHM (or very 
similar) equations showing that at large time the flow is strongly organized 
and dominated by structures [H], [12], [IB], [H] (and references therein). At 
the oposite limit the turbulent regime can be treated with renormalization 
group methods [15]. 
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The formulation of the problem exclusively in terms of experimentally- 
accessible quantities, velocity and vorticity seems to not allow too much free- 
dom in elaborating a theoretical model from which the stationary states to be 
determined. In all our considerations we will be guided by the analogous ex- 
perience in the case of the Euler fluid. In that case the existence of a parallel 
formulation (although the mathematical equivalence is still not fully proven) 
has been decissive. That model consists of the discrete set of point-like vor- 
tices evolving in plane due to a potential given as the natural logarithm of the 
relative distance. Comparing with the differential equation of the Euler fluid 
this alternative formulation provides something fundamentally new: it splits 
the dynamics into two objects of distinct nature: point-like vortices and po- 
tential, or, in other words, matter and field carrying the interaction. Going 
along this model we are led to consider the standard treatment in these cases. 
We must assume that each of these objects can evolve freely {i.e. indepen- 
dent of the other) and in addition there is an interaction between them. This 
formulation is standard in electrodynamics (classical or quantum). Going to 
continuum, the discrete set of point-like vortices becomes a field of "matter" 
(which must be assumed in general complex), the potential will be a free field 
(called "gauge" , similar with the free electromagnetic field) and, in addition, 
there is the interaction between the matter and the gauge fields (similar to 
the classical jfiA^^ charge-field interaction). Basic properties of the original, 
i.e. physical, model impose constraints to this field-theoretical formulation. 
The presence of vorticity requires a particular form of the potential between 
discrete vortices: it is the curl of a sum of natural logarithms. This can only 
be derived from a gauge-field Lagrangean density of Chern-Simons type, in- 
stead of Maxwell type. Therefore the field-theoretical model will contain the 
Chern-Simons Lagrangean. The matter field has a nonlinear self-interaction 
that reflects the stationary structure of the free matter field. The coupling 
of the two fields (the interaction) is minimal, via the covariant derivatives. 
In a previous work we have formulated this model and have obtained, in this 
way, a purely analytic derivation of the sm/i-Poisson equation for the Euler 
fluid [IS]. 

This paper will develop a model for the Charney-Hasegawa-Mima equa- 
tion along the same lines. However several differences will impose new fea- 
tures of the model. 

In a previous work [17] we have investigated a model based on the simi- 
larity with the superfiuid field theory, the Abelian-Higgs model. This is able 
to describe positive fluid vortices. We have also proposed, without details, a 
more extended model, which seemed able to describe the physical vortices of 
plasma and atmosphere within the regime of the CHM equation. 
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The full development of this model is the main objective of the present 
work. 



2 The physical problem and the Charney-Hasegawa- 
Mima equation 

The analytical model we develop in the present work is intended to describe 
a system characterized by the following elements: 

1. the existence of an intrinsic length; this is the Larmor radius ps for the 
two-dimensional plasma immersed in a strong, confining, transversal 
magnetic field; and the Rossby radius for the two-dimensional quasi- 
geostrophical approximation of the planetary atmosphere; 

2. the existence of a condensate of vorticity as a background in the system, 
in the absence of any perturbation. This background consists of the 
cyclotronic gyration of ions, for the plasma; and of the Coriolis effect 
resulting from the planetary rotation, for the atmosphere. 

These basic elements are very general and other systems belong to the 
class that is defined by them. In particular the non-neutral plasma produced 
in laboratory experiments and the vortices produced on a fluid in a rotating 
tank (see Schecter, Nezlin, Hopfinger and Van Heijst, etc.) 

To make the discussion more specific we will refer in the following to the 
physical model developed for the two-dimensional plasma and atmosphere, 
which leads to the equation of Charney-Hasegawa-Mima (CHM): 




(90 dip 
dt dy 



[(-Vx0 X n) ■ V±] Vl0 = 



(1) 



where ney = — n x V^lnnQ. For simplicity we will refer to the problem of 
plasma physics, the adaptation to the problem in the physics of atmosphere 
being easily done (see [E]). The quantities appearing in the Eq.([T]) are the 
physical ones after having been normalized : 



{x,y) 




(2) 



t = t^^y'n, 
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where Qci = l^lB/rrii, pg = Cg/Qci, = Te/rrii. The derivation of the 
equation, in the drift instabihty in tokamaks is done in the Appendix A. 
For comparison, the Euler equation is 



~dt 



or (3) 



^ (vi0) + [(-Vx0 X fi) • V J vl0 = 

where (j) is the streamfunction. 

The similarity between the two equation, Eq.([T]) and Eq.(I3]) is apparent. 
However, with regard to the form and properties of the stationary states 
of the two equations we should not attempt of simply taking any previous 
conclusion derived from the Euler fluid context into the CHM context. This 
is because the naive stationarity {d/dt = 0) imposed in the two equations 
leads to an equation with a vast degree of generality and does not provide, 
by itself, a clear identification of the final vortex shapes. We may represent 
the family of all possible solutions of the naive stationary limit as a subset 
in a function space. In the time evolution the two equations produces two 
series of configurations representing functions belonging to two distinct paths, 
ending in this set at distinct points {i.e. configurations). 

Compared with the equation for the ideal fluid (Euler equation) the CHM 
eq. is not scale invariant p^. To see this we make the rescaling of the space 
variables 

{x,y) {x',y') = {\x,Xy) 
Expressing the Euler equation in the new variables, we have 

I {Vl<f>) + X' [{-V^<j> X n) . V'J Vl<j^ = 

The factor can be absorbed in a rescaling of the time variable and the 
equation preserves its form. By contrast the equation CHM becomes (for 
simplicity we take k, = 0) 

(1 - X'Vl) ^-X' [(-V'^0 X n) . V'J V20 = 

While the factor can be absorbed by time rescaling, the factor A^ in the 
first paranthesis cannot be absorbed. The form of the equation is invariant 
only for A = 1 which means the space is measured in units of pg. The equation 
CHM exhibits an intrinsic spatial scale, which is Ps- 
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3 An equivalent discrete model 



In the case of the Euler fluid, there is an equivalent model whose dynamical 
evolution is considered to be identical with that of the physical system (for 
a list of references, see [16]). It consists of a collection of discrete point-like 
vortices in two- dimensions evolving from mutual interaction defined in terms 
of a potential. The potential is the natural logarithm of the relative distance 
of two vortices. This model has been proposed and used long ago (Kirchoff, 
Onsager, etc.) but the rigurous proof of the equivalence between it and the 
physical Euler description is a difficult mathematical problem [19J. 

For the CHM equation there is a similar model: a collection of point-like 
vortices interacting by a potential that has a short range. This model has 
been proposed in meteorology by Morikawa [20] and Stewart [21] (see Horton 
and Hasegawa [H]). For a set of N vortices with strength Uj, j = 1,N 
with instantaneous positions the streamf unction ip {x, y) has the following 
expression 

^ (r) = ^ (r) = J2 ^jKo (m |r - I) (4) 
j j 

where m is a constant. The differential equation from which the contributions 
to the streamfunction ip in Eq.(jlj) are derived, is 

(A - m^) ijj (r) = -2'KUJj6 (r - r^) (5) 

in two dimensions. For a single vortex of strength u placed in the origin, the 
azimuthal component of the velocity can be derived from the streamfunction 

V0 = — = —umKi (mr) (6) 
or 

At small distances ^ 

f 61 ~ - for r ^ (7) 
r 

which is the same as for the Euler case, where ip (r) is given by the natural 
logarithm. The streamfunction decays fast at large r since the modified 
Bessel function of the second kind Kq decays exponentially at large argument 

tp ~ —;= exp (—mr) for r —>■ oo (8) 
\/mr 



This means that the vortices are shielded. The elementary vortex of the 
Charney-Hasegawa-Mima equation is localised by m~^ and one can associate 
a finite spatial extension to it, ps in plasma physics, pg in the physics of 
atmosphere [14] . 
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The equations of motion for the vortex ujk at {xk,yk) under the effect of 
the others are 120j 




dW 



(9) 



dyk 
dW 

dxk 



where 



N N 



(10) 



i=i j=i 



This is the Kirchhoff function for the system of interacting point-hke vortices 
in plane. It is the Hamihonian for the system of vortices. If we introduce 
the versor of the normal to the plane, n, the equations can be expressed 



We will develop a formalism for this discrete system which in the contin- 
uum becomes a field theory. 

4 The main components and the main steps 
of construction of the continuum model 

We will develop a continuum model whose equations of motion could repro- 
duce, in the discrete approximation, the Eqs.(jl]), © and the energy ffTUl) . The 
model must be defined in terms of a Lagrangean density for two interacting 
fields: 

• the field associated with the density of point-like vortices (f){x,y); we 
will call it the matter (or scalar, or Higgs) field; and 

• the field associated with the potential carrying the interaction between 
the vortices; we will call it the gauge field. 

4.1 The gauge field 

We note that the interaction potential Kq (|r — rj|) (assume that m is nor- 
malised as m = 1) appearing in the discrete model proposed in meteorology 
is similar to the potential appearing in the Euler problem, ln(|r — rj|) in 



(11) 
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the following sense: they both have topological properties, in sense to be 
explained. For the Euler equation, the potential can be represented using 
the angle made by the line connecting the reference {rj) and the current (r) 
points with a fixed line, and in order to remove the multivaluedness one has 
to make a cut in the plane from the center (where is singular) to infinity 



Euler fluid : ^ = (-V0 x nf (12) 

(Here j label the point-like vortices and a , (3 label the coordinates of the 
position vectors, rf, a = 1,2). Since 

e-/^— = e"/3^^1nr (13) 



we see that the potential in Eg. (1121) is expressed through the Green function 
of the 2D Laplacian, defined by the equation 

V^lnr = 27r5^(r) (14) 

The potential is obtained by applying the rotational operator e'^^dp on this 
Green function. 

The CHM case is similar, with the difference that the Kq {mr) is the 
Green function of the Helmholtz operator, as results from Eq.([n]). The series 
representation of the function Kq is 

oo 2k 

^oW = -/o(.)ln| + X:^^^(A: + l) 

we see that close to the origin the two potentials are similar 

Ko{r^O) = -ln^ + ... (15) 

We note that the potential in the Euler fluid case may be presented as a 
singular pure gauge 



—5-^^ = - — Aarctan^ (16) 
2n 277 ar" x 



2ti dr 
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The individual contributions to the potential are the derivatives of the angle 
9 made by the particle position vector with an arbitrary fixed direction in 
plane. In the case of the CHM equation we have from Eqs.(!ll), ([9]), (IHl) 

CHM plasma : — ^ = (-V0 x ii)" (17) 

(Jjb 

= e°''^y^ujj ^ [m\r — Tjl Ki{m\r — Tjl)] 

~^ |r — r,-| 

Using the small argument expansion (see formula 8.446 in [23] or, alterna- 
tively, the Eq.([IS])) 



m 



I r — r j I i^i (m I r — I ) — > 1 for | r — r j | ^ 



we note that Eq.( |T7|) can be written like Eq.( |T6l) . The function in the right 
paranthesis only changes the spatial decay. 

The angle 6 {x, y) is a (multivalued) scalar function and we have in both 
cases a typical situation of the type 

(-V0 X nf ~ dj (18) 

so that the potential can be considered at large distances a pure gauge 

(-V0 X n)" ~ g-^dg (19) 

with g E U (1) i.e. g = exp {i9). In other words, for every point {x,y) on a 
large circle on the plane, we have a value of the angle 6. These potentials 
have therefore a topological nature, since they map the circle at infinity in 2D, 
(r oo) onto the set of values of the angle 6, which is also a circle. This is 
a typical homotopic classification of states and the potentials are classified 
into distinct sets characterised by an integer, representing how many times 
the circle in the plane is covered by the circle representing the values of 6. 

The main difference between the Euler and the CHM cases is the short 
range of the potential in the latter case. If we use formally the concept of 
photon (a "particle" that mediates the gauge interaction as in electrodynam- 
ics), one can say that in the Euler case we have the usual (two dimensional) 
massless photon, whereas in the CHM case we have a massive photon. The 
fact that the photon is massive is another way to express the fast spatial de- 
cay of the potential function, i.e. the short range and we will often use this 
formulation. From physical reasons we know that this short spacial range 
must be of the order of ps, the intrinsic length in the CHM equation. The 
need that the equations of motion lead to a short-range potential (finite-mass 
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photon) and that the potential has a topological nature represent constraints 
for the part of the Lagrangean density coming from the gauge field. 

As we have shown in the case of Euler fluid [16], the topological nature of 
the potential is provided by the Chern-Simons (CS) Lagrangean 

Ccs = '^e^^^A^df.A^ (20) 

where e"'^'^ is the totally antisymmetric tensor in 2 + 1 dimensions (a, (3 and 
7 can take three values: 0, 1, 2, corresponding to the time and the two coor- 
dinates (x, y)) and k is a constant. This Lagrangean is essentially the density 
of "magnetic" helicity. It is known that this Lagrangean does not lead by 
itself to dynamical equations for the potential since it is first order in the 
time derivatives; it only represents a constraint on the dynamics, analogous 
to the Lorentz force in an external magnetic field given by the combination 
of K with the other constants of the model. We can heuristically say that the 
Chern-Simons Lagrangean induce vortical effects on the dynamics. This will 
become more clear later. 

The gauge field dynamics can be introduced either by coupling the Chern- 
Simons potential with the matter field, or by including the Maxwell La- 
grangean density 

Cm = -If.^F^" (21) 

or both. We note that any of these combinations provide a finite mass for 
the photon but the way they do that is different. 

For the combination Maxwell and Chern-Simons, the vortical effect of the 
Chern-Simons part induces decay of the field on structures of small spatial 
scale, with an extension governed by the "external magnetic field", k. This 
corresponds to the gyration motion. 

For the combination Maxwell and matter field the generation of the pho- 
ton mass is due to a classical Higgs mechanism. The matter field has a 
nonlinear self-interaction which vanishes at certain non-zero values of the 
field. Therefore the extremum of the action implies the minimum of this 
self-interaction potential and the matter field will take one of these values 
(called: the vacuum value) at infinity. This is the symmetry breaking leading 
to the Higgs mechanism. The motion of the photon in a polarisable medium 
consisting of this background matter density induces a finite mass effect for 
the photon. Then the value of the mass for the photon (the short range of 
the spatial decay) is determined by this vacuum value of the matter field and 
by the coefficient of the Maxwell contribution in the Lagrangean (the electric 
charge). 
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The combination Maxwell, Chern-Simons and matter has therefore two 
possible ways to obtain a finite mass for the photon: the gyration (due to 
Chern-Simons) and the Higgs mechanism, due to the finite background of 
matter field corresponding to one of its "vacuum" values. There is a mixing 
of these two ways and there are two possible masses, or short ranges for the 
gauge potential. At the limit where the Maxwell term is suppressed from the 
Lagrangean, the short range of Chern-Simons with matter is recovered. 

We argue that including the Maxwell term is not necessary in our model 
describing the flow governed by the CHM equation. This would only provide 
for the gauge field an independent dynamical evolution since, even when the 
matter field is absent, the Maxwell Lagrangean leads to plane waves, i.e. a 
propagating field, without any meaning or justification in our case. 

4.2 The matter field 

The matter field (f) is associated to (without being identical with) the density 
of point-like vortices in the discrete model. The matter field must be com- 
plex since the vorticity carried by any point- like vortex appears as a sort of 
electrical charge (only complex fields can represent charged particles). The 
kinematical part of the matter field in the Lagrangean consists as usual in the 
squared momentum but with the covariant derivatives, to refiect the so-called 
minimal coupling with the gauge field 

Ckin=-l{Dy)\D^<P) (22) 

where 

D^ = d^ + (23) 

In the Hamiltonian formulation of the discrete vortices model for the 
Euler equation it has been derived an equation connecting the gauge field 
with the "density" of the point-like vortices. Going to the continuum {i.e. 
field-theoretical) version it appeared that the only way to keep this constraint 
was to assume a self-interaction of the matter field. The same reasons act 
in our present case, but now the problem is more complicated. The self- 
interaction potential V (0) must have a minimum at a nonzero value of the 
matter field such as to ensure the background that will induce (together with 
the CS term) the short-range of the gauge field. Comparing with the Euler 
case (we neglect the various constant factors) 

VEuleri<l>)-^\<P\'' (24) 

the simplest form would be 

VcHM{(l>)-^{\(l>f-v^Y (25) 
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where v is the vacuum value of the matter field. However, it will be shown 
below that this form cannot provide, for the Lagrangean density consisting of 
Chern-Simons and matter part, the most symmetrical extremum of the action 
functional for the system. This particular state is called self-duality and we 
adopt the point of view that this is a fundamental requirement on the model. 
In particular in our previous paper [16j it was shown that the stationary 
states of the ideal fluid obeying the sm/i-Poisson equation correspond to the 
self-duality. 

In a different context [26], [27] the form of the self-interaction potential 
able to support self-duality has been found as 



and we will work with this form. 

4.3 The necessity to consider "pairing" of the fields 

The main physical content of the states governed by the CHM equation is the 
vorticity, which is organizing in large scale vortical structures. The vorticity 
and the velocity fields generate the kinematic helicity, which is a topological 
invariant for a dissipationless fiuid/plasma. It has been shown [22] that the 
helicity is determined from the boundary condition (this becomes evident 
in the Clebsch representation of the velocity): the values of helicity on the 
boundary of the volume are sufficient to determine the value at any internal 
point. We understand that besides the fields defined in the internal points, 
we need to consider fields that carry the information from the boundary 
toward the interior, on equal foot with fields that carry informations from 
the points of the internal volume to the boundary. This pairing of functions 
suggests that all quantities involved in our model will be matrices. The model 
becomes non-Abelian and the quantities are elements of the algebra of the 
group SU (2). 

A more formal explanation of the necessity to adopt a non-Abelian alge- 
braic structure of the theory results from the consideration of the spinorial 
nature of the elementary point-like vortices and from the Parity, Charge 
Conjugation and Time inversion invariances of the theory. 

5 The field theoretical formalism 

The continuum limit of the system of discrete point-like vortices is a field 
theory. From the discussion of the previous Section, we have the field the- 




(26) 
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oretical model: covariant, SU (2), Chern-Simons for the gauge field and 6*^ 
order self-interaction for the matter field. 

• gauge field, with "potential" A^, (/i = 0, 1,2 for {t,x,y)) described by 
the Chern-Simons Lagrangean; 

• matter ("Higgs" or "scalar") field described by the covariant kine- 
matic Lagrangean {i.e. covariant derivatives, implementing the mini- 
mal coupling of the gauge and matter fields) 

• matter-field self-interaction given by a potential V (0, 0^) with 6*'' 
power of 0; 

• the matter and gauge fields belong to the adjoint representation of the 
algebra SU (2) 

The Lagrangean density for such a model has been used in (2 + 1) field 
theories and reads 

C = -KE^-'PtT (^d.A.Ap+'^A.A^A,^ (27) 



-tr 



(D^0)t (D^0) 



-r(0,0t) 

What follows is already exposed in field-theoretical literature, in particular 
in Dunne [30], [35], [31]. For the Abelian version, see [28] . 

The transformations of the space and time variables must be connected 
through the condition of the general covariance of the theory. The metric of 
the space-time is 

/-I 0\ 

9>.u = 9^''= \ 10 (28) 

V 1/ 

This means that we have to take account of the covariant and contravariant 
coordinates of vectors, tensors and operators. We can use both notations 
{Ai,A2) and {A^,Ay) since no confusion is possible. We have 

= it,x,y) (29) 

Af" = {A°,A\A^) (30) 
A^ = {Ao,A,,A2)=g,uA''={-A'',A\A') 
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the derivation operator is 



(31) 



= g'^'d^ (32) 
d d d 



and 



dt dx dy ^ 
The covariant derivatives are 

D, = d^ + [A,,] (33) 

(note that we need not introduce an electric charge, e). We write the detailed 
expression 

D^cj) = g'"'D^(j) (34) 
For comparison we write them in detail (see also Eq. (l30l) ) 

= + ^00 - Mo, 1^ + Ai0 - (j)Ai, ^+A2(f)- (j)A2 

[ot ox ay 



y ot ox Oy J 



(35) 



The Hermitean conjugate of a matrix is the transpose matrix with complex 
conjugated entries. For Eq.i lMj) the Hermitian conjugate is 

(D»^ = (9»^ + [A^,0]^ (36) 



or, in detail 



{Dyy = |-^ + 0t^ot_^ot^t^ (37) 
ox 



+ _ + 0tA2t_A2t0t} 



dy 
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keeping the following rules (* is complex conjugate and T is transpose) 



This means 



{A 



(3^ 



*T 



(39) 



A; = 1,2 



It has been found that the only possiblity this model has to reach self- 
dual states is to choose a matter field nonlinear self-interaction given by a 
sixth order potential [26] 



V 



) 



4^2 



:tr 



([[0,0t]^0] [0,0^,0] -^V) . (40) 



The trace is taken in a finite dimensional representation of the compact 
simple Lie algebra Q to which the gauge field A^ and the charged matter 
field and 0^ belong. 

The Euler Lagrange equations are 



KB 



tJ.p 



dV 

90t 



where the field strength is 

Fii„ = d^Ay — dyA^ + \A^, Ay 



p,u 



(41) 
(42) 

(43) 



(We may note that the second equation of motion shows a typical Chern- 
Simons' action property, i. e. the proportionality between the magnetic field 
and the current density, which in physical systems would be called force- 
free). 

The current is 



(^^ 



(44) 



with the conservation (covariant) 



D^'J,, = 



(45) 



17 



The Gauss law constraint is the 0-component of the second equation of 
motion 

_ ^ (^012^^2 + ^021^21) = iJ' (46) 



—iJi 







since 



J' = 9'' J, = 9'' Jo = -Jo (47) 
Using the Eg. fl44l) we have the Gauss law 

2«:Fi2= [<p\Do<P] - [{Do<P)\<p\ (48) 

in the nonabelian form. The equation written above is the 0-th component 
1^^012 _ ^021 _ q£ ^Yie equation of motions connecting the field tensor 
with the current of matter, Eq. fl42p . We identify in the right hand side the 
magnetic field, since 

Ff,, = d^A,-d,A^ + [A^,A,] (49) 






-Ey 


-Ex 


Ey 





-B 


Ex 


B 






(/i, 1/ = 0, 1, 2) i.e. 

F12 = -B (50) 
The energy density of the system is 

£ = tr ((Z}o0)U/^o0)) + (51) 
+tr f (Z^fc0)^ {Dtq 



+V (0,0t) 

The energy density can be rewritten in the Bogomolnyi form, i.e. as a 
sum of squares plus a quantity that integrated over the plane becomes a 
lower bound for the energy. The space part of the Lagrangean containing 
covariant derivatives can be written [37], [35] 



tr ((Dfc0)t {Dk<p)) = tr ((/^-0)^ (I)-0)) - ^tr (0^ [F12, 0]) (52) 
with the notation 

D± = Di± 1D2 (53) 
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In the following we will verify the equality 



- ^tr (0t [Fi2, 



2k 



-tr 



[[0,0t] ,0]^Do0- ,0] Po0)^) (54) 



Using Eq. (H8i) we replace F12 and obtain 
-ztr (0t [Fi2,0]) 
= —iti 



(55) 



^([0t,Do0] - [(/^o0)^0]),0 
^tr{0t|[0t^/)^0]0_ [(Do0)^0]0 



tr 

2k 



t,Do0l+0^(/^o0)^0 



2k 



tr{0t [<p\Do<P\<P 

(^o0)^0 



+ 



(^o0)^0]} 



We expand the commutators in order to collect together the factors of Dqc, 
and respectively (Z)o0)^- 



-ztr [F] 



12, 



(56) 



2k 



tr{0t (0t (Z}o0) _ (Z}o0) 0t)0 



-0t ((Do0)^0- 0(^00)^)0 
+0Vf(Do0)t0_0(Do0)t 
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-rtr(0t[Fi2,0]) (57) 

(L»o0)^00 

+0^(^00)^0 

-(/)t#t(Do0) 

+(/.V(Do</')</'^ 

+(/)V(L>o</')^0 
-(/)t#(L>o0)^} 

We take separately the terms containing Dq4> and use the cychc symmetry 
of the Trace operator 

(/.Vt(Do0)0 (58) 
-0t (Do^ 



0V)}(Do0) 
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and analogously for the factors of (Dq^)^- 

-</.t (Do0) (59) 

+</.V(^o0)V 

+0V(^o0)V 
-0t00(Do0)t 

-0V0(Do0)t 

= {(00t_0t0)0_0(00t_0t0)jp^0)t 



It results 

-^tr(0t[Fl2,0]) (60) 
= -^tr{[[0,0t] ,0t] (Do0)+ ,0] (Do0)^} 

and we will prove that 

tr{-[[</.,0t]^0t]}=tr{[[0,0t]^0]^} (61) 
The right hand side is 

(62) 



and the left hand side 



-[[0,0t],0t] (63) 
-{(00t_0t0) 0t_0t (00t_0t0)} 
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and one can see the identity of the two expressions. Then 



-ztr (0t [Fi2,0]) 

^tr { [ [0, 0t] , 0] t _ [ [0, 0t] , 0] 1 1 



(64) 



This is the expression that is used in Eq. fl52|) 



We have 



tr (^{Dk(py iDk(p)) 



tr - ztr (0t [Fi2, 0]) (65) 

tr ((Z}_0)^ (D_</.)) 

+ ^tr I [ [0, 0t] ^ 0] t (^^0) _ [ 0t] ^ 0] p^0)t I 



6 The term containing the time-derivatives 

We remind that the objective of this calculation is to find the configurations 
for which the energy Eq. flSlj) is minimum. According to the usual approach 
to this problem, we will try to reexpress Eg. (1511) as a sum of squared terms 
plus a "residual" term. The squared contributions, being always positive, are 
minimum when the corresponding expressions are zero, while the additional 
term in general has a topological meaning. In the particular case of the CHM 
fluids, the residual energy cannot be associated to a topological quantity, for 
reasons that will be discussed later. Since we are no more guided by the 
physical significance of the additional energy as resulting from a topological 
property of the system we must accept that there is no unique way of sepa- 
rating in Eq. (l51l) the square terms and the additional term. We present in 
the following two such formulations and discuss them comparatively. 

6.1 First mode of separating the squared terms in the 
energy expression 

Now we have to write in detail the term from the Lagrangian containing 
the zero-th covariant derivatives. This is done by including the constant v 
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representing the asymptotic value of the charged field. 

tr f (Do0)^ {Do<p)) (66) 

t 

2k 



tr| (Do0-^([[0,0t],0]-^;20) 



^ ... fi\\M M^ ..2 At 



\tr(([[0,0t],0]_^2^)t([[^,^t],^]_^2^) 



4^2 

This expression is obtained after several cicliv translations of the factors, 
which is allowed under the Trace operator. The last term, when we return 
to Eq. llST]) . cancels the potential V (0, 0''") given in (HOl) . We now understand 
that only this choice allows to write the energy in the Bogomolnyi form. 

The total energy in the system, Eq. flSlI) . written in Bogomolnyi form, 
results from the expressions ( l65l) and ( !66l) 

£ = tr|^(^Z}o0-^([[0,0T,0]-t^V)y (67) 



+tr ((/^-</))^ (D-<j 
+ ^tr (0t(Do0)-(Z^o0)V 



The energy contains a sum of positive quantities (squares) and is minimised 
by those states where these terms are vanishing. The last term shows that 
there is a lower bound the energy 



• 2 

£ > ^tr (2)^0) _ (/)^0)t 0j (68) 

6.1.1 The first form of the self-duality equations 

The vanishing of the squared terms in the energy leads to the self-duality 
equations 

D_<p = (69) 
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Combining these two equations such as to put in evidence the gauge field 
whose expression is in Eq. dC.lTl) . (the calculation is presented in detail 
in Appendix C) 

= (70) 

rh 

Using Eqs.f l69l) we can derive a new expression for the energy in the self-dual 
state 

Sso = f^tr{<l>^{v'<P -[[<!> A^] A])) (71) 
which is the saturated lower bound shown above. 

6.2 Second mode of separating squared terms in the 
expression of the energy 

We look for an expression of a square term that differs from the previous one 
by a change of sign within the first two terms of Eq. fl66|) 

tr I (^Do<P + ^ ( [ [0, , 0] - ^ V) y (^/^o0 + ^ ( [ [0, 0T , 0] - ^^0) ) I (72) 

= tr |((Do0)^ - ^ ([[0,0t] ,0] -^V)') (^00+ ^ ([[0,0T ,0] -t;^))} 
= tr|(Z}o0)^(/^o0) 

([[0,0^,0] -^^'0)^(^00) 
+ (Do0)^f ([[0,0^,0] -^V) 
+ 4^([[0,0^0]--V)^([[0,0T,0]--^0)} 
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The two median lines are expanded and the full expression is rewritten 
tr I (/^o0 + ^ ( [ [0, 01 , <f\ - ^ V) ) ' (/^o0 + ^ ( [ [0, 01 , 0] - ^ V) ) I 

( [ [0>T )^ (^00) + ^^'0^ (^00) 
+ ^ (Do0)^ ( [ [0, 01 ' 0] ) - ^ (^00)^ 

+^ ( [ [0> 01 > 0] - ^'0) U [ [0, 01 , 0] - ^^0) } (73) 
and we can now get an expression for the first contribution to the energy 
tr|(Do0)U^o0)} 

= tr I (^Do0 + ^ ( [ [0, 0l , 0] - V) y (^00 + ^ ( [ [0, 0l , 0] - ^^V) ) 
-tr I ( [ [0, 0l , 0] ) ^ (1^00) + y/^^ (^00) 
+ ^ W)^ ([[0,01 >]) - ^ (^00)^ ^V} 

-tr{4^([[0>1>]--^0r([[0>1>0]--^0)} 

With the two expressions detailed above, we can rewrite the energy 
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Eq.dSH) 

8 = tr ((Do0)U^o0)) +tr (75) 



-tr ([[0,0^] ,0])^ Po0) + ^t;20^ W) 

+^ W)M [ [0> 0T > 0] ) - ^ 

-tr{i^([[0,0^0]--V)^([[0,0T,0]-^^^0) 
+tr (Z}_0)) + ^tr I [ [0, 0t] , 0] ^ (Do0) - [ [0, 0T , 0] (/^o0)^ 

We note the cancellation of the term which is the potential V (0^, 0) and we 
obtain 



£ = tr|(^Do0 + ^([[0,0T,0]-^V)J (^/^o0 + ^([[0,0T,0l 
+tr 

(76) 

where the additional term in the energy expression is 

= -trl^-^{[[<P,<P^,<p]yiDo<P) + ^v'<PHDo<P) (77) 

+^(Do0)U[0>0T,0] -^Po0)^'0 



-tr 



{[[0,0t],0]^(Do0)-[ [0,0^,0] (1^00)^} 

{[[0,0t],0]^(Z}o0)} 
— tr |(Do0)^ [[0; 0^] ; 0] | (uslng cyclic permutation in Trace) 
|-tr|0t (Do0) - (Do0)V} 
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We write 

ga ^ ga(l) ^ ^a{2) ^ ^a(3) j^yg^j 

for the last three hnes of the equation above and these will be calculated 
below. At this point it is more useful to focus on the set of equations at 
self-duality that are derived from this choice adopted in Eq. (17^ . 



6.2.1 Second form of the Self-Duality equations 

The NEW equations as they result from the alternative Bogomolny form of 
the action are 







(79) 



^00 = -—{[[ct>A^]A]-v^) 
We introduce the field tensor and it is shown in Appendix C that 



^0) 



- \ -I 
K 



[0t,z}o0] - [(I^o0)^0]) 



Now, we use the NEW equations 



(A 



^1) 



2k 



Inserting the two operators from the equations above, and after finding that 
the two commutators in the Eq.f lC.6l) are equal and opposite, we get an 
expression for F+_ as 



m 



where we replace the NEW expression of DQ(j) , i.e. Eq.( l8T|) . obtaining 



2i 



(83) 



-4 [<P\[[<PA']A]-v'<p] 
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Then the NEW equations at self-duahty are 







(84) 



1 



+ 



We note that this form of the self-duahty equations differs from the pre- 
vious one by the opposite sign of the right-hand side term of the second 
equation. 

Using the definition of Eq. (|53l) the left hand side of the first equation of 
motion (HTll can be written 



7 The group theoretical ansatz 

7.1 Elements of the SU (2) algebra structure 

In order to solve the self-duality equations it is considered, as in the case of 
the Euler equation, the Lie algebra of the group SU (2). Then the Chevalley 
basis is |32| 




(85) 




tr {E+E_) 



2 



H 



±2E, 
1 



(86) 




E± are step (ladder) operators 



The 2x2 representation is 




(87) 




(88) 




(89) 
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The Hemitian conjugates of the generators are the transposed complex con- 
jugated matrices 

El = E_ (90) 
eI = E+ 
= H 

These adjoint generators will be used to express the adjoint fields in the 
calculations based on a particular ansatz. 

7.2 The fields within the algebraic ansatz 

According to Dunne [31], [33], the following ansatz can be adopted 



(j) = Y^(l)aEa + (t)-ME.M (91) 
a=l 

since the rank of SU (2) is r = 1. We take the Hermitian conjugate, which is 

0t = 0*4 + 0^^1 (92) 
= <PlE^ + (j);E+ 

In this ansatz the matter Higgs field is represented by a linear combination of 
the ladder generators plus the generator associated with minus the maximal 
root. 

The gauge potential is taken as 

A+ = aH (93) 
A_ = -a*H 

The notations with + and — correspond to the combinations of the x and 
y components, with the coefficient i for the y component. 

7.2.1 The explicit form of the equations with the ansatz 

The gauge field tensor 

F+_ = d+A_-d^A+ + [A+,A_] (94) 
= d+ {-a*H) - 9_ (aH) + [aH, -a*H] 
= i-d^a* -d^a)H 
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where 

d+ = d, + idy^2-^ (95) 

d 

d- = dx- idy = 2— 
oz 

We will have to calculate, with this ansatz, the terms of the equations. 
The right hand side of the second equation 

(96) 

will be calculated using the commutator 

^ [cP,E+ + <P2E.,riE- + <P;E+] (97) 
= Ml[E+,E_]+Ml[E_,E_] 
[E+,E+] + <P2<p*2 [E-,E+] 

[0,0t] = ^l^,[E^,E_] (98) 
+0;02 [E-.E+] 

= {4>14>, - 4>14>2) H (99) 
= (Pi - P2) ^ 
where we have introduced the notations 

Pi = 101 1' (100) 

P2 = |02P 

The next step is to calculate 

[[4>, 0^] , 0] = [(pi - P2) H, ct>iE+ + ct>2EJ\ (101) 

This is 

[[0>t],0] = (pi-p2)0l[/^,£^+] + (pl-p2)02[/^,£^-] (102) 
= 2 (pi - P2) (01^+ - 02^-) 

The next level in the commutator is 

,0] = + 1;'02^- (103) 

-2 (pi - P2) + 2 (pi - P2) 02£;- 
= PE+ + QE_ 
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where 



P = ^2^^_2(pi-p2)0i (104) 
Q = ^7^2 + 2 (pi - Pa) 02 

Returning to the Eq. (l70p . the full right hand side term is 

[v'^ - [[0, 0t] , 4>] , 0t] = [PE+ + QE_, 4>IE_ + 4>;E^] (105) 

= P<Pl[E+,E_] + Q^;[E_,E^] 
= {P<Pl-Q<p;)[E+,E^] 

= {p<pi-Q<p;)H 



or 



= { {v^<Pi - 2 (pi - P2) 0i) 0t - {V^<P2 + 2 (pi - P2) 02) H 

= (pi - P2)) pi - + 2 (pi - P2)) P2} H 

7.3 Using the algebraic ansatz in the first version of 
the SD equations 

The second self-duality equation Eq. (!70l) becomes, using Eqs.lIMl) and fll06p 

-£-|^4('"-''^)^'^-2(P.+P.)] (107) 
Now we turn to the first self-duality equation 

D_0 = (108) 
and its adjoint form. It has been defined 

D_=Di- iD2 (109) 

then 

D_0 = 1^ + [A,, Ct>\-l^-l [Ay, 0] (110) 

To proceed further we express the components of the potential 

= A^ + iAy = aH (111) 
y4_ = A^- iAy = -a*H 
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Then 



l{a-a*)H (112) 



Ay = l-{a + a*)H 



Then 



D (j) ^ -i^) E++ (— -i—] E (113) 

+l(a-a*)0i [H,E+] 
+l(a-a*)02 [H,E_] 
-i^(a + a*)0i 
_il(a + a*)(/.2 

L>_0 = (^^-z^ + 2i(a-a*)0i-2i(a + a*)0i^£;+ (114) 

+ (,9^"^9^"^2^'^-'^^'^^ + 22^'^ + '^^'^V 
= 

From the exphcit form of the ladder generators we obtain the equations 
derived from the first self-duahty equation 

|i-i|^ + 20.a- = O (116) 



7.3.1 The explicit form of the adjoint equations with the algebraic 
ansatz 

Now we consider the adjoint equation (also derived from the extremum of 
the corresponding part of the action expressed in the Bogomolnyi form) 

(D_0)^ = (117) 
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We have 



D 



dx 



(118) 



where the adjoint is taken for any matrix as the transpose complex conju- 
gated. The change of the order in the commutators is due to the property 
that for any two matrices Ri and R2 the Hermitian conjugate of their com- 
mutator is 

[i?i,i?2]^ = {RiR2-R2Ri)^ (119) 
= {R2R1 - R1R2) 

(* is complex conjugate and ^ is the transpose operators) and we take into 
account that in the expression of 0^ we have already used the Hermitian 
conjugated matrices of E±. 

The Hermitian conjugates of the gauge field matrices are 



Then 

Dl 

We recall that 
The we have 



4 = \{a*-a)H^^\{a*-a)H 

Al = -l(a* + a)i7t^_J_(«* + „)^ 



| + ,^+l(«*_«)[,^]_l(«* + «)[,^] 



(t>'^^(t>\E_ + (t>*2E+ 
a) IH] 



— + i— + -{a 
dx dy 2 



X 



dx 



dy 



E_ + 



dx 



+ i 



1 
2 

dy 



a* + a) [, H] 



E^ 



+^{a*-a)<p; [E+,H] 
-^{a* + a)<Pl [E_,H] 



(120) 



(121) 
(122) 

(123) 
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or 

= 2^E, + 2^E^ (124) 

V g^* dz* ^ ' 

+^(a*-a)0t (2E_) 
+1 (a* - a) 0; i-2E^) 
-^{a* + a)<Pl (2E_) 
_l(a* + a)0;(-2E+) 

The equation becomes 

= (2^ + {a*-a)<j>l-{a* + a)<f>l]E. (125) 



9^* 



dz* 



+ { 2^ - (a* - a) 0; + (a* + a) 0^ ) 



Here we have made use of the identifications 

d ■ 9 _ ^ c 
dx dy dz 



d d d , , 

^ + ^^ = 2— (126) 



and 

d . d ^ d ^2 

dx dy dz 

The resulting equations are 

d6* 

2--i - 2a</)t = (128) 
dz* ^ 

and 

2|^ + 2a0; = O (129) 
which represent the adjoints of the first set, Eqsf lllSI) . as expected. 

7.3.2 Using the two sets of equations 

Now we consider the first equations {i.e. those refering to 0i) in the two sets, 
Eqs.([II5]) and ffT^ 

2^-2a*0i = (130) 
dz 

2^-2a0i = 
34 



From here we obtain the expressions of a and a* 

a=Ain(0*) (131) 

a* = ^ In (00 (132) 
The left hand side of the second self-duahty equation fll07p is 

da* da d d , , ^ d d , , , , 

= -2a^|l"to) + ln(«l 

dzdz* ^^^^^ ) 
In the differential operator we recognize the Laplacean, 



A 



dzdz* 



Equating the expressions that we have obtained for the left hand side and 
respectively for right hand side of the second self-duality equation (11071) we 
obtain 

- i A Inpi = (pi - p2) [2 (pi + p2) - v^] (134) 

The second equations (those refering to 02) in the two sets Eqs. (11161) and 
(11291) give the result 

a* = In 02 (135) 
dz 

and 

a = -^ln0* (136) 

from where we obtain the form of the right hand side in the second self-duality 
equation, (I107p 

= 2g^|ln (,#.,) + In (,^;)] 
dzdz* ^ 
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The final form is 



:i38i 



Tfie riglit liand side in Eqs. fll34p and (11381) is tlie same and if we substract 
tfie equations we obtain 








(139) 



A In pi + A In p2 
A In (P1P2) 

Tlie function In (P1P2) is an arbitrary harmonic function and this aspect will 
be discussed later. For the moment we simply take a constant, convenient 
for normalization, 

P1P2 = v^/lQ (140) 
With this relation we return to the equation for pi, (I134p 

ih = I n^ I z I /;i -H 



Pi 



2 Pi + 



Pi 



:i4r 



We add the zero- valued Laplacean of a constant to the left side and factorise 
in the right hand side 



lAlnpi-lAln(t;V4)=4^ 



Now we introduce a single variable 



Pi 



f2/4 



Pi 



Pi I 

Pi 



t;2/4 



P = 



Pi 
t;2/4 



P2 



and obtain 

1a, 1 
-Alnp = - — 

2 ^ 4 V 
We make the substitution 



i/j = \np 



and we obtain 



1 

^ ( — 1 [exp(^) -exp{-tp)] 



X < - [exp {ip) + exp (-^Z")] - 1 
1 fv^Y 

- ( — 1 sinh ip (cosh — 1) 



(142) 
(143) 

(144) 
(145) 

(146) 
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(^-^j A?/" - sinh?/^ (cosh?/' - 1) = (147) 

Exactly the same equation would have been obtained starting from the one 
for p2, (11381) after a change of the unknown function, ip —ip. 
After normalizing the coordinates by the length /t/f^, we obtain 

AV^ - sinh?/' (cosh?/; - 1) = (148) 

This is the equation governing the stationary states of the CHM equation, 
resulting from the first form of the SD equations. 



7.3.3 Calculation of the additional energy for the first version of 
the SD equations 

We start from the energy as integral of the density of the Hamiltonian 
Eq.f lSlI) . Since all other terms in the expression of the energy are posi- 
tive (they vanish after adopting the self-duality and the particular order 
potential), the energy is bounded from below 



IV 



S^^—tT{(j)HDo(j))-{Do^ 



where the second of the equations at self-duality Eq. 



IS 



(149) 



(150) 



and 



Then we have 



2k 



Ah? 



tr{0t ([[0,0t] ^cf\-v^) + ([[0,0t] ,0] (151) 



and we will prove that the second term in the curly brackets is equal with 
the first. 



t 



(152) 



0t [0,0t]0_ [0,0t] 0t0_^2^t0 
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We can apply in the second term (underlined), the cyclic symmetry of the tr 
operator, moving successively the factors and 0^ in the first position and 
have 

tr {([[</), 0t],0] -v'<py<p] (153) 

= tr{0^ [0,0^] [0,0^ 
= tr{0t[[0,0t],0] -t;V0} 
= tr{0t([[0,0t],0] -v'<P)} 

and the equality with the first term in Eg. (115 II) is proved. It results 

S ^ -i^tr {0t ([[0,0t] ,0] -t;^)} (154) 

but at self-duality (since we have already used the equations derived from 
self-duality) the limit is saturated 

Ssn = ^tT{<l>^v'<p- [[0,0T,0])} (155) 

We can obtain the explicit formula using the algebraic representation of 
the fields 

= 0iE+ + 02E_ (156) 

0t = <pie. + <p;e+ 

and recall the previous result 

v^(P - [[0, 0t] , 0] = PE+ + QE^ (157) 

where 

P = „2^^_2(pi-p2)01 (158) 
Q = 17^02 + 2 (pi - P2) 02 

A detailed calculation, starting from Eq. (ll55p : 

£sD = ;J^tr(0t(t;20- [[0,0t],0])) (159) 
= ^tTmE^ + <l);E+)iPE+ + QE^)} 
= — tr{0tP^-^+ + <PIQE^E^ + <PIPE+E+ + (PIQE+E^} 
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We calculate separately 
tr {(I)IPE_E+) 



= tr [0^ (^;Vi-2(pi-p2)0i)^-^+] (160) 
= Vi - 2 (pi - P2) Pi] tr 1^ |j J j 

= v^pi - 2 (pi - P2) Pi 

= tr [0t (t;V2 + 2 (pi - P2) 02) ^-^-] (161) 

= 0t (^2^2 + 2 (pi - P2) 02) tr ° ° 

= 



tr(0;p£;+£;+) = 



tr (^;^0i - 2 (pi - P2) 0i) ^+^+] 

02 (^^01 - 2 (Pl - P2) ( Q Q 



(162) 



tr(0;QE+E_) 



= tr (r;202 + 2 (pi - P2) 02) E+E_] 
= [^^2 + 2 (pi - P2) P2] tr J Q 
= ^;^P2 + 2 (pi - P2) p2 



Summing up the contributions 



^SD — 



2^2 



2^2 



Expressed in the normalised variable, the energy is 

,.2 

^SD 



V 

2«;2 



[v^ (pi + P2) - 2 (pi - P2) ] 

^2 1 



4p 



- 2 P-r - 



2 2 1 \ 2 



2«:2 V 4 



11/ 1 

p 2 V P 



4 4 p 

2" 



1 / 1 
p 

p 



(163) 



[v'^pi - 2 (pi - P2) Pl + v'^p2 + 2 (pi - P2) P2] (164) 
[v^ (pi + P2) - 2 (pi - P2)^] 



(165) 
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Introducing the streamfunction p = exp (ip) we get 



S, 



SD 



\2 (sinh ijj) — 2 cosh 'ip] 

8 Pi 

— — — [(cosh?/') — cosh-?/; — l] 

Ps 



(166) 



or 



SsD = v"^-^-: \- (cosh?/))^ + cosh?/) + ll 

pi^ 



(167) 



We note that this expression must be integrated over the plane (the factor 
1/pI will ensure the correct dimension) and the dimension of the energy is 
actually given by = Qci- 



7.4 Using the algebraic ansatz in the second version of 
the SD equations 

We now turn to the second version of the SD equations ( l84l) and introduce 
the algebraic ansatz. Then the second equation of the second version of the 
Self-Duality becomes 

F+- = (168) 
= -\[v^-2{p, + p2)]{pi-p2)H 



Using Eq.(l94D 

- 1^ - 1^ = (Pi - P2) K - 2 (Pi + P2)] (169) 

OX+ OX- K"^ 

From the first equation of self duality, which is common to the two choices 

a=Ain(0*) (170) 

a* = ^ In (00 (171) 
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The left hand side of the second self-duahty equation fll07p is 



da* 



da 



dz* dz 



2 — —\ri. 

dz* dz 

d^ 



2——\n\ 

dz dz* 



-2 
-2 



dzdz* 
d^ 

dzdz* 



[ln(0i)+ln(0t)] 
ln(|0i|') 



Since the Laplace operator is defined as 



A = 4 



d^ 



dzdz* 



we get 



;Aln(|0i|^) 



(pi - P2) [v'^ - 2 (pi + P2 



or 



^Alnpi = ^(pi-p2) [v'^-2{pi + p2 



and we can now replace 



oAlnp = - 



2\ 2 



P2 



4 



1 

P - - 
P 



Alnp 



,4 /I 



P 



K 

V 

2 2^ 
and introduce the streamfunction 

p = exp (ifj) 

,2X 2 



4-2 p + - 

pyj 



1 



p 



or 



V 
2\ 2 



sinh ip {1 — cosh ■?/') 



Alp + I — ) sinh?/) (cosh?/' — 1) = 



The unit of space is 



1 

Ps 



V 
K 



(172) 



(173) 

(174) 
(175) 

(176) 
(177) 



(178) 
(179) 

(180) 
(181) 



and the equation results 

Alp + sinh ij (cosh ^ ~ 1) = (182) 

All the other calculations, in particular those implying the function 02 
and the complex conjugated, (pl and (p2 are similar to the calculations made 
for the first version of the SD equations. 
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7.4.1 Calculation of the additional energy for the second version 
of the self-duality 

The additional term in the Bogomolnyi form of the energy, in the second 
version, Eq. ( 1781) consists of three contributions. The first contribution is 

£'^(i)^ltr{[[0,0t],0]t(^^^)| (183) 

and we use the previously derived expression 

[[0, 0t] , 0] = 2 (pi - p2) {ct>,E+ - <f)2E.) (184) 

and 

[ [0, 0^ , 0] ^ = 2 (pi - p2) {<PIE_ - <P;E+) (185) 
Also we use the following relation 

v^^-[[(l),^^],^]=PE+ + QE^ (186) 

where 

P = ^2^^ _ 2 (pi - p2) 01 (187) 
Q = v'^(f)2 + 2 (pi - Pa) 02 

Using the second (new) equation of self-duality we have 
gai^) ^ ltr{[[0,0t],0]tp^0)| (188) 

= ^tr 1 2 (pi - p2) {<PIE^ - r,E+) (-^ ( [ [0, 0T , 0] - v'<p) 



^2 (pi - p2) tr mE. - <P;E+) {-PE^ -QE.)} 



The trace is 



tr {(0tE_ - 0;E+) {-PE^ - QE_)} (189) 

0t (-P) tr {E^E+) (trace is 1) 

+ (-02) i-P) (^+P+) (trace is 0) 

+0t (-Q) tr (E^E^) (trace is 0) 

+ (-02) i-Q) tr (P+P-) (trace is 1) 

0t (-P) + (-0;) i-Q) 

-0t [t;20i - 2 (pi - P2) 0i] + 0; [t;^02 + 2 (pi - P2) 02] 
-pi [v"^ - 2 (pi - P2)] + P2 [v^ + 2 (pi - P2) 
-t;^ (pi - P2) + 2 (pi - P2) (pi + P2) 
- (pi - P2) [^^^ - 2 (pi + P2) 
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and the first contribution to the residual energy becomes 

S''^'^ = l(pi-p2){-(pi-P2)b'-2(pi+p2)]} (190) 



Now we calculate the second contribution to the residual energy 

£<'^ = -^tr{(Do</')U['/'>T>]} (191) 
= -^tr |(Do0)^ 2 (pi - P2) (01^+ - 02^-)} 
and we replace, according to the new second SD equation 

{Do^y = (^-^^ ([[0,0t],0]_^2^)t (192) 

= ^(-p^+-gE_)t 

= -A(p*£;_ + Q*£;^) 
and replacing in the previous equation 

S"^'^ = -^tr{(L'o(/))^2(pi-p2)(0i£;+-(/'2£^-)} (193) 
= -^tr I (^-^ (P*E_ + g*E+)^ 2 (pi - P2) (0ii^+ - hE-) 
^ ;2 (pi - P2) tr {(P*£;_ + Q*E+) {(f)^E+ - (f)2E_)} 



2k^ 

The trace is calculated separately 

tr{(P*E_ + Q*E+) {(t)iE+ - (t)2E^)} (194) 
= P*0itr{E_E+} (trace is 1) 

+Q*0itr {E+E+} (trace is 0) 

+P* (-02) tr (trace is 0) 

+Q* (-02) tr {E+E_] (trace is 1) 
= P*0i + Q*(-02) 

= [wVl - 2 (pi - P2) 0l] * 01 - [w^02 + 2 (pi - P2) 02] * 02 

= pi [v^ - 2 (pi - P2)] - p2 y + 2 (pi - P2)] 
= -y^ (pi - P2) - 2 (pi - P2) (pi + P2) 
= (pi - P2) [-y^ - 2 (pi + P2)] 
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and the second contribution is 

S""^'^ = --,ipi-p2){pi-p2)[v'-2ip, + p,)] (195) 

= (pi -P2)^ K -2(pi + P2)] 

We note that the two contributions S"'^^^ and f are equal. 
The third contribution is 

^"^'^ = -^tr {<P^ (^00) - (i^o0)^ 0} (196) 
We can use the new second SD equation 

Do<l>^-^{[[cl>,cl>^],cl>]-v'cl>) (197) 



together with 
We find, as in previous cases, 



v^(l) - [ [0, 0t] , 0] = + QE_ (198) 



and 



i^o0 = -^{-PE+-QE^) (199) 

= f (P£;+ + Q£;_) 



(L>o0)^ = ^{PE+ + QE_)^ (200) 
= ^(P*E_ + Q*E+) 
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Then 

S"^'^ = -|^tr{0t(L»o0)- W)^^} (201) 



~tr l^t^ + QE_) - ^ {P^E_ + Q*E^) cj> 

iv^ ( i 



= — {0*Ptr (P_^+) (trace is 1) 

+0iQtr (E_E_) (trace is 0) 
+02^tr (£^+£'+) (trace is 0) 
+0;Qtr {E^E_) (trace is 1) 
+P*0itr [E_E^) (trace is 1) 
+Q*0itr [E^E^] (trace is 0) 
+P*02tr [E_E_) (trace is 0) 
+g*02tr(E+E_) (trace is 1) 

We obtain 

^"^'^ = -J^{0tp + p*0i + 0;g + gV2} (202) 



2^2 {<^t - 2 (Pi - P2) 0i] + 0; [^^'02 + 2 (pi - P2) 02] } 
— {pi y - 2 (pi - P2)] + P2 y + 2 (pi - P2)] } 

[v^ (pi + P2) - 2 (pi - P2)^] 



2^2 

£:a(3) ^ i_ [^2 + _ 2 (p, - p^f\ (203) 



2«;2 

Now we collect all results 



ga ^ £a(2) ^ ^a(3) ^204) 



(Pi - P2)^ - 2 (pi + P2)] 
(Pi - P2)^ - 2 (pi + P2)] 

^2 

+ ^[^^' (Pl + P2)-2(P1-P2)'] 
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there are three powers of v and we separate the coefficients 

v° : (205) 
^ (Pi - P2f 2 (pi + P2) + ^ (Pi - P2f 2 (pi + P2) 

= ^ (Pl - P2)^ (Pl + P2) 



(206) 



P2)' 



(Pl - P2)^ 



P2)' 



P2)' 



: (207) 

^(Pl+P2) 



and the total expression is 



4 3^^ V 

^"-^ (P^ - P'^f (P^ + ^^2) - ^ (Pi - P2)' + (Pi + P2) (208) 

= ^ (/'I - Z'^)' [4 (pi + P2) - 3v'] + ^ (Pi + P2) (209) 
Introducing the normahzation 

P = ~V77 = (210) 

^2/4 P2 




V 



16^2 

,6 



P- 



P 



V" / 1 



P + 



-3 



(212) 
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V 



,6 



■a 




2 ^ 

(sinh ip) [2 cosh ip — 3] + -—r cosh 



(213) 



As in the case of the Eg. (11671) we will note that the expression must be 
integrated over the plane, which removes from the coefficient the dimensional 
factor 



8 Discussion on the versions of the SD equa- 



The possibility of formulating the expression of the energy as a sum of squared 
terms plus a (additional) term with topological significance (known as Bogo- 
molnyi formulation) is fundamental for the self-duality. In our case, the CHM 
fluid/plasma cannot associate a topological significance to the additional en- 
ergy, a characteristic signalized by Lee 1991 and by Dunne. This induces a 
certain imprecision in the choice of the way of separating the squared terms, 
with consequences on the form of the equations, etc. This aspect will only 
be discussed briefiy here, with the only intention to compare few possible 
choices. 

We have shown how to derive two versions of writting the total energy 
of the system as a sum of squared terms plus an additional {residual) term, 
while this one has no topological significance. After adopting the algebraic 
ansatz we arrive at two different equations for the scalar function ip which 
we associate with the physical streamfunction of the CHM fluid. 

The choice of the version that has the correct physical signiflcance should 
be done on the basis of the supersymmetric invariance of the extended fleld 
theoretical model. However, even from this advanced point of view, we can 
expect at most an indication which will not be applicable directly to our 
problem. The fleld theroetical model for the CHM equation shows signiflcant 
differences compared with topological theories in 2D. The SD equations, in 




and the dimension of the energy is given by = Q, 



tions 
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both versions, lead to time dependent solutions, therefore the stationarity 
tj^ical for the solutions obtained from the Bogomonlyi form in other theories 
is here lost. The topological aspect is also lost, the additional energy, in both 
versions, is not proportional with the total winding number induced by the 
vortices present in the plane. We must note however that the point of view 
that results from the SUSY extension of the theory (a natural extension) 
favorizes the Eg. (11481) or possibly, as mentioned below, the Abelian version, 
Eq. fl214p . This even if the additional energy does not have a topological 
meaning. 

One possible help comes from looking at the theory of the CHM fluid as 
being a development of the theory for the Euler fluid. Since the asymptotic 
states of the latter are goverend by the sm/i-Poisson equation, we expect that 
the nonlinear term of the equation for CHM fluid to be of a similar nature. 
Or, we see that the first set of SD equations leads to a sign which is oppo- 
site to the one appearing in the sm/i-Poisson equation. The physical form 
and properties of the solutions would be, in that case, completely different. 
Or, one expects that an ideal fluid without an intrinsic length will transform 
smootly into a fluid which has an intrinsic length, as the CHM fluid. Al- 
though this is not an argument, we take this as a sort of indication that the 
second choice is more appropriate and we adopt Eq. fll82p as the equation 
governing the asymptotic states of the CHM fluid. 



Integrand of E (1/4)[cosh(i|/) - (cosh{\|/))^ +1] 
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Figure 1: The integrand of the energy Esd which is the additional term in 
the first Bogomolnyi form. 

A different approach can be developed on the basis of the analysis made 
by Lee (1991) of the first set of equations at SD. After adopting the algebraic 
ansatz, Lee finds that the boundary conditions on the two functions (^\ and 
02 lead to restrictive choices : basically it results that the second ladder 
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Integrand of E (1/4) [2(cosh(\[i))-' - 3(cosh(v/)) -(cosh(\|f)) +3] 




ignitude of the streamfunction \|/ 



Figure 2: The integrand of the energy £sd which is the additional term in 
the second Bogomolnyi form. 



generator should not be present in the algebraic form of 0. Then the equation 
resulting from the first set of SD equations is identical with the equation 
which is derived from the Abelian version of the theory 



Alp = exp {2iIj) — exp (ip) 



(214) 



Contrary to the Equation 01481) . Eq. (l214p has physically interesting solutions, 
consisting of rings of vorticity. The way this theory can arise from the Abelian 
dominance in the full field-theoretical (FT) model of the CHM fluid and the 
physical consequences of this theory will be analysed elsewhere. 

Another possible criterion for the choice of one or another Bogomolnyi 
form is the behavior of the additional energy. In some cases, the solution 
of the equations at self-duality {i.e. by taking to zero the squared terms) 
show a higher energy due to the residual term, when compared with other 
choices. In particular we have found that the following form has apparently 
lower energy for many of the solutions of its associated self-duality equation, 
in practical applications like the tropical cyclone, etc.: 



Alp + - sinh ip (cosh ip — 1) = 



(215) 



or a version in which the harmonic function is a constant different of 1, i.e. 
taking 



P 



Pi 



(4p) 



instead of Eq. 01761) we have 



Alp H sinh ip (cosh ip — p) = 

2p 



(216) 



(217) 
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This form arises by taking 



in the Eq.( l75i) . with 



A 



and has an associated residual, non-topological energy 



1 

X - 

4 



11 3 

— {smhip) (— 2 + cosh-?/;) H — cosh?/' 



which is shown in the Fig. 



Integrand of E (1/4) [(1 1/8)sinn(v)^(-2+cosh(\|f)+(3/8)cosh(v)] 




Figure 3: The integrand of the energy Ssd which is the additional term in 
the third (for A = 1/2) Bogomolnyi form. 



9 Various forms of the equation 



The crucial step in the derivation of the form Eq.f ll48p or Eq.f ll82p is the 
choice of a solution for the equation (11391) 



Aln(pip2) =0 



(218) 



with the consequence that the product pip2 is the exponential of a harmonic 
function. This leads to the conclusion that the asymptotic states of the CHM 
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fluid can be described by a class of differential equations, parametrized by 
the harmonic functions. Without developing this aspect here we mention 
that the fact that instead one equation we find a class of equations has a 
simple physical significance. The difference resides in the background mo- 
tion and the equations from this class describes vortices on a background of 
fluid/plasma motion which has zero physical vorticity. Then the "stream- 
functions" ipi and ■02 that can be introduced separately for pi respectively 
for p2 will differ by a function whose Laplacean is zero, i.e. a potential flow. 
We note that by physical vorticity we understand the vorticity perturbation 
above the condensate background, i.e. above Qci- 

In the following we mention few examples illustrating this freedom of 
choice of the zero- vorticity background flow. 

The separation of variables in polar coordinates 

A$ = (219) 

^(r,e) = R(r)e(e) 

dPR 1 dR 



with 



+ r— -Z^^ = (220) 
+ n^Q = 



dr'^ r dr 



d9^ 



9.1 Solution 1 

The first choice was 

with the consequence 



ln(pip2) = (221) 
Pip2 = 1 (222) 



9.2 Solution 2 

A different (almost arbitrary) choice 

A/i = (223) 

h = exp (x) sin (y) (224) 

Then 

In (P1P2) = exp (x) sin (y) (225) 
P1P2 = exp [exp {x) sin {y)] (226) 
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9.3 Solution 3 (general solution in cylindrical coordi 
nates) 

A different cfioice for tlie cylindrical harmonic function $ = In {pip2) 

1 



X [C cos {nO) + D sin {nO)] 



and for n = 0, 



$(r,^) = {AO + B) 

X [Chir + D] 



(227) 



(228) 



9.3.1 Example polar 1 

For example, 

$ (r, = ar cos + 6 
Consider the choice with a particular value for h 

In (P1P2) = cos 6* + 6 
P1P2 = exp (or cos ^ + 6) 



(229) 

(230) 
(231) 



- exp I ar cos 1 



16p' 

for p a positive constant. Define p as before 

Pi 



P = 



l^exp [\ar cos 9) 
1 (\ 



— — exp -ar cos 9 
P2 4p ^ V 2 



-A In 



2 
1 



p— exp -ar cosd' 



= ^ (Pi - P2) [2 (pi + P2) - v"^] 

AC 

where 

= -^(pi-p2) [2(pi+p2)-^'] 



1 

X I2 



■;r / 1 
p— exp I -ar cos ( 
4p 



exp I -ar cos ( 

p 4p ^ V 2 



/I ,\ it;^ /I , 
p — exp -ar cos exp -arcosU 



(232) 



(233) 



(234) 
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RHS 



— exp -ar cos 



P-- 
P 



(235) 



V 

X — exp 
4p 



or cos ^ I < 2 ( p + 



4pexp 



-ar cos ^ 



1 fv 



2\ 2 



\ n 

X (p - - 
V P 



— ) exp (ar cos 9) 



(236) 



X 



1 ^ M 

- \ P -\ — \ — P exp — ar cos o' 



and the Right Hand Side 
LHS = 



-A In 



p— exp -ar cos 
% ^V2 



(237) 



1 fv^ 
— Alnp + Aln — 

2 V4p 



1 A n 

— A -ar cos/ 
2 V2 



The equation is 

1 



-Ap 



(238) 



4p2 \ n 



,,2\ 2 



exp (ar cos ^) I p 



pexp 



However, since we will note 



-ar cos 9 



(239) 



it would have been easier to do that from the beginning. 
The equation 



exp (i/j) = ^ 



Pi 



exp (iar cos^) 



(240) 



4p 

1 (\ 



— exp ( -ar cos 9 

P2 4p V 2 
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generates explicit forms for pi^2 that will be inserted in the equat 

Pi — — exp I ip + -ar cos v 



4p " V 2 



with the same result. 

The equation now reads 

— A'^ H — - exp (ar cos 6*) sinh 

p2 



where the space unit is now k/v"^ = Ps- 

9.3.2 Example polar 2 

The second simple polar choice of A$ = 

$ (r, 6*) = a In r + 6 

ln(pip2) = alnr + 6 

P1P2 = exp (a In r + b) 

= exp (b) r" 

The normalization suggets 



_ Pi 



4p' 

P2^P 



p4p 
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The equation is 



jAlnp-^Aln(r«/2) 



1 

xl2 



4p pip 



— V 



(248) 



-Alnp 



(249) 



1 fv 



,2 \ 2 



4p^ V / V P 



X 



2^P + -J-pr 



-a/2 



After introducing the substitution Inp = ip and measuring the space in units 
k/v'^ — ps we have 



or, better 



A-^ H — -r^sinh-^ fcosh-^ — — ^ ~ 



— A'^ H — - sinh ^0 (r" cosh ^0 — p) = 



(250) 



(251) 



10 Discussion on the physical meaning of the 
model 



10.1 The short range of the potential 

It is considered that the scalar field is very close to the vacuum value 

We calculate the current in the region of vanishing space- variation. 

= -i {0^ {dy + [A^ 0]) - (a> + [A^ 0]) 0t 

- (9^0t + [0t, A^t] ) + (5^0t + [0t^ ^Mt] ) I 
= -2 {0t (9^0) - (^^^0) 0t _ (5^0t) + (5^0t) 

+0t [A^ 0] - 0] 0t - [0t, ^^t] + [0t^ ^Mt] I 



(252) 
(253) 
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Since we consider that the field is almost constant (and equal to v) we can 
negelct all terms on the first line and obtain 



Let us consider the explicit expressions for the fields 



A. 



A,, — 



A^ 



l{a-a*)H , = ha*-a)H 
^(a + a*)H , A^^ = (a* + a) H 



and 



(f) = (f)iEj^ + (t)2E- 
Then, for x, the first part of the formula for the current is 



[a — a 



(a - a*) [(t)lE_ + - 02^-] 

(a - a*) {0^01 - 0^02 E_]} 

(a -a*) (-|0i|'/J-|( 
- (a - a*) (pi + P2) 



and the second part 

[0, [0^^i] 



0t£;_ + 0;E+,-(a*-a)i/ 



= -(a*-a)[0,0ni?-,i/] + 0;[i?+,i/]] 

= ^(a*-a)[0i£;+ + 02£;_,0t2£;_-0;2£;, 
= (a* - a) {0i0t[^+,^_]- 020; [£;_,£;+]} 

= (a*-a)(|0i|'//+|02|'//) 
= {a* -a){px + P2)H 

The X component of the current is 



7 ~ 



{- (a - a*) (pi + P2) + {a* - a) (pi + P2) H] 
2i{a - a*) (pi + P2) H 



(254) 



(255) 



(256) 



(257) 



(258) 



(259) 
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For the far regions we take the value 

Pi + P2 ~ 'y^/S 
We return to the potential notation, (a — a*)H — 2A^ 

Analogous for the y component of the current, 



and the second part 
[0, [0t,^.t]] 



it, 1 (a + a*) (01 [H,E+]+<p, [H,E_]) 



= -ia + a*)[(f>\i<f>^2E+-h'^E_)] 

= -i(a + a*) + -02^-] 

^ _i (a + a*) {0^01 E+] - 0^02 ^-]} 

= -i(a + a*) (-|0i|'i/- |02|'i/) 

= i (a + a*) (pi + P2) H 



^lE_ + ct>lE+,--{a* + a)H 



-- {a* + a) [0, <Pl [E_, H] + 0^ [E+, H]] 



1 

Yi 



{a* + a) [0i^+ + 02E_,0t2^_ 



= t{a* + a) {0i0t E_] - 020; 
= ^(a* + a) (|0i|'/J+ |02|'i/) 
= i (a* + a) (pi + P2) -f^ 

The X component of the current is 

Jy ~ -i {i {a + a*) {pi + P2) H + i {a* + a) {pi + P2) H} 
= 2{a + a*)ipi + p2)H 

As before, we replace here pi + p2 ~ "^^ and {a + a*)H = 2iAy, 

We take the temporal component of the potential in the form 

A° = bH 



0\*T 



(260) 
(261) 

(262) 



(263) 



(264) 



(265) 
(266) 



b*H 
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Then the temporal component of the current density is (cf. Eg. (12541) ) 

J°^-z{[0t,[A°,0]] + [0,[0t,^ot]]| (267) 

and we calculate again the terms, with the particular choice Eg. f l266p 

= [(l>\h{(l>,[H,E+] + (^2[H,EJ^)\ (268) 

= 6 (0i2E+ - 022E_)] 

= 26[0*E_ + -02^_] 

= 2h{-\ct>r\^H-\(l>2fH) 
= -2b{p, + p2)H 

and the second part 

[0, [0t,^ot]] ^ + (269) 

= b*[<P,<Pl[E^,H]+<p;[E+,H]] 
= b* + 02^-, 0t2E_ - (P;2E+] 

= 2b* - Ml [E^, E+]} 

= 2b*{\(PifH+\(l)2fH) 
= 2b*ipi+p2)H 

The result for the current in the approximation of the constant background 
is 

J° ^ -^{[<P\[A',<P]] + [<P,[<P\A'^]]} (270) 
= -i{-2b{pi+p2)H + 2b*{pi+p2)H} 
= 2i{b-b*){p^ + p2)H 

We then obtain 

J° = 2iv^A'' (271) 

Then 

J'^ = (J°, J", J^) = 2iv'' [A\ A\ Ay] (272) 

or 

J'^ ~ 2it;2A^ (273) 

With this value of the current density we return to eguation connecting 
the gauge field with the matter current. The eguation is 

- KE^^'PF^p = iJ" (274) 
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To transfer the antisymmetric tensor e'^'^^ in the other side, we multiply by 
e^o-T and sum over repeated indices 



'^fiar'^ ^ up 



up 



K 



(275) 



err Tcr 



K 

i 

2k' 



-f 



A direct relation with the previous result is obtained taking the explicit form 
of from Eg. (12731) 



^ err r\ ^ UfJT <^ 

2k 

Introducing the expression for the field 

Far = dijA^ — d^-Acr 

we apply the derivative operator dr and sum over the index r 

= F) TP- 

„ liar '-'t "J 

2k, 
i 



-—e^^rdr{2iv''A'') (from Eq. (12731)) 



The term on the right hand side is 
and here we replace, from Eg. (I276p 



2^aTpJ- T 



9^ ^rariA^ 



(276) 



(277) 



(278) 



(279) 



(280) 
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Further the product of the two antisymmetric tensors e is expanded 

drdaAr - drdrA„ = S f.^rOrA'' {{lom Eq.^Wj) (281) 

1 

= -T^^ar^K (from Eq.(l279D) 

= -T^^ar^g^'^-SravA'^ (from Eq.(^) 

1 / 2\ 2 

1 / f ^ 



2 ( ^ ) 9^ ^uT^^TarjA^ 



The exphcit expression for the sum 



9^ ^uT^^rar] (282) 

9^^ (^o-lO^lOTy + £^cr20^20»7) 
~^9^^ i^aOl^Olri + £^o-2l£^21r?) 
+9^^ (^(T02^02r? + ^erU^Ur]) 

(5a25r,2 + 5,l5^,l) 
+6'^^ ('5a2'5»72 + 5f70'^»?o) 

+^(T2'^r?2 + 5a0^r)Q + ^al^-ql + '^o-0<^??0 



— 25o-O^J70 

This is replaced in the Eg. (12811) 



1 / 2\ 2 



drd^Aj- — drdrAfj — 2 ( "/7 ' 9^°'^<7Tfi^TarjA'^ 



1 

2 V K 



2\ 2 



-) 5.nA° 

2\ 2 



K ■ 



— 1 ^ctqAq 



Taking the gauge condition 



drA^ = (283) 
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it results the equation 

drdrAo -[-] Ao = (284) 



2N 2 



The solution of this equation is, in cylindrical geometry, 

(r) = Kq [mr) 
From here we conclude that the mass of the photon is 



m = — (285) 

K 

and this mass is generated via the Higgs mechanism adapted to the Chern- 
Simons action. The photon gets a mass because it moves in a background 
where the scalar field is equal with the vacuum value, nonzero value. 
From physical considerations, we know that 

1 

m = — = — (286) 
10.2 A bound on the energy 

We can express the total energy of the system as the space integral of the 
time-time component of the energy- momentum tensor 



Stot = J d^T''' (287) 
A useful formula (Gradshtein 6.561 formulal6, [23]) is 



l^^x^K. (ax) dx = 2-^a--r (^^f^) T (^^f^ 



(2^ 



for 



Re{p + l±u) > (289) 
Rea > 



Then, for = 1 and a = 1, u = 0, 

POO 

/ xKo{x)dx = [T{l)f 
Jo 



(290) 
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(cf. Gradshtein 8.338). This must be used with Eq. fnUj) to calculate the total 
energy of a system of vortices in plane. 



W 



cont 



277 j Sroj'^K {m\r I — r2\) 

POO 

/ (mr) d [mr) Kq [mr] 



(291) 



, UJ 



m 



Then we have that the 2D integral over the plane of the continuum version 
of the energy of a system with discrete vortices is constant multiplying the 
square of the elementary quantity of vorticity, which was before associated to 
each elementary vortex. This corresponds actually to the value of the energy 
in the field theoretical model, precisely at the self-dual limit, Eq.(!68l). 



10.3 Calculation of the flux of the "magnetic field" 
through the plane 

Start with the second differential equation of self-duality Eq. flTOj) 

1 



This has been calculated previously, with the result 

yct> - [[0, , 0] , = {v'-2 (pi + P2)) (Pl - P2) H 

Then 

F+- = -^{v^-2 (pi + p2)) (pi - p2) H 
We return to Eq. lHSl) (the Gauss law constraint) 



We can express in detail this constraint, using the Eqs.(??) and (??) 



(292) 

(293) 
(294) 

(295) 



Ft 2 — 

2k 



r,E. + <p;E^,~—{PE+ + QE^] 
Zk 



(296) 



-—(PE+ + QE_) ) ,<P,E+ + <P2E- 
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The first term is 



<l>lE_ + <l>*E+,-—{PE+ + QE_) 



I 

'2k 



{<l>lP[E_,E+] + <f>lQ[E+,E^]) 



2k 



{cf>lP-cf>lQ)H 



— {v^-2{p,+p2)) {Pi-P2)H 



The second term is 



Then 



-— {PE+ + QE^) ) , (j>^E+ + (t>2E. 



= —[P*E_ + Q*E+,(j>iE+ + (t>2EA 
Zk 



2k 



= -—{P*ct>i-Q*<t>2)H 

= -^(^'-2(P1+P2))(pi-P2)i^ 



i^i2 = —i^—{v^-2{p, + p-,)){p,-p^)H 
+^ - 2 (pi + P2)) (Pi - P2) 



The result gives us the magnetic field 



F12 — —B 



Comparing with 



2^2 



(t;2-2(pi+P2)) {Pi-P2)H 



= -{v^-2{p^ + P2)){pi-p2)H 



4^ 



P 



P + 



P 



H 
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we note the relation 



Fi2 -2^+ 



The flux is 



d\ (^;2-2(pl+p2)) (P1-P2) 



(302) 



(303) 



The quantities pi and p2 are not normalized, therefore it is preferable to 
change to the variable p 

(304) 



We note that 



v'^/A P2 

tr (00^) = pi + p2 

it] = {p^-p2)H 



(305) 



tr 



Pi+P2= {vVa) p + ^ = -{p+- 



(pi - P2) if 



£(,4,. 



The flux is normalised as 



d\ [v^ - 2 (pi + P2)) (pi - P2) 
2 



v'-2lp— + —- 



V 



II f ,2 f 1 



d'rlp-- 



4 4p/JV4 4p 
1 / 1 



2l^^P 



P— - — - 
1 



(306) 



(307) 



NOTE. In the abelian relativistic model the following relation exists 
between the flux and the minimum energy [30] 



S 



Abelian 
SD 



— tr 
2 



which means 



r Abelian 



2 + 



(30^ 
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END OF THE NOTE 

We calculate the scalar field self-interaction potential 

2 At 



^ tr ([[0>t]^0] ([[0,0^,0] -^V) (309) 



tr 



{PE+ + QE_)^ {PE+ + QE_) 



1 

^tr [{p*E. + Q*£;+) {PE+ + g£;_)] 

^ [P*Ptr + g*gtr 



{p*p + Q*g) 



since the other traces are zero. Using the notations pi and p2 we have 
1 



V 



1 
1 



{P*P + Q*Q) (310) 
2 - 2 (pi - p2)] ' pi + [^;' + 2 (pi - p2)] ' P2} 

= 4^ ^'"^P^ ~ ^'"'^ ~ ^2) Pi + 4 (pi - P2)^ Pi 
+t'V2 + 41;^ (pi - P2) P2 + 4 (pi - p2)^ P2} 

^ {y^ (pi + P2) - 4i;^ (pi - p2f + 4 (pi - p2)^ (pi + P2)} 



We can express the potential in terms of the normalised variable 

V (0, (/.t) 

{f^^ (Pi + P2) - (pi - P2)^ + 4 (pi - p2)^ (pi + P2) } 



(311) 



4^2 

1 

4«;2 



f ^ f ^ 1 

^T + Tp 



t.2 1 



-4v^ p 

^' A A p 

4 4 p 



+4 p 



V V 1 

^T + Tp 
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V{4>,4>^) (312) 




P2 p, 

Figure 4: The potential 1^(0,0^). The points where pip2 = 1 are shown. 

'jm') 




Figure 5: Same as figure H] with a different view. 

10.4 Comment on the possible associations between 
the field-theoretical variables and physical vari- 
ables 

The field theoretical model has been developed as the continuum version of 
the system of discrete vortices interacting via a short range potantial. On 
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the other hand, the original model, represented by the CHM equation, is 
expressed in terms of three variables 



V = —V X e^^ip 
OJ = V X V = Alp 

(where is the versor perpendicular on the plane). These variables have a 
clear physical meaning. We would like to understand the possible connection 
between the variables of the field theory (j^, A^, F^^, cp) and parameters 
[k and v"^) and the physical variables. We make the observation that the 
conservation law for the scalar and vorticity field [H] is (including explicitely 
the normalisation factor ps) 

jj d^r [i)'^ + [p^ipf] = const 

We make an integration by parts of the second term 

d'r {Vtpy = 1 1 d^r (V^ ■ V^/-) = 

rfV [V {ip ■ Vip) - ijAi)] 

The integration of the first term is transformed 

d'^rV {ip -Vip) = (pd\-{ipVij) 



Since we have the definition v = —Vip x we see that the scalar product 
excludes the component of Vtp that is normal to the contour, which is a 
circle of very large radius. The diamagnetic effect makes that the velocity 
of gyration of particles on Larmor orbits generates a macroscopic velocity 
of the fluid which is tangent to the circle, therefore V^p has only a nonzero 
component, normal to the circle. We conclude that there is no contribution 
from the first term in the integrand. Then we have 

J J d'^rip (^tp — p'^A'ip) = const 

In particular this means that a vacuum state, of zero energy, corresponds to 

plA^P -^p = 
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or 

p2^-^ = (313) 

On the other hand, we have a characterisation of the vacuum state in the 
field theoretical model, obtained as the asymptotic state of the fields at large 
distance. There the scalar field is almost constant and the space derivatives 
are vanishing (this is also known as the large wavelength approximation). It 
will be shown below that the current and the potential verify the 
relation 

f - 2iv^Af' ~ (314) 

If these two relations Eqs. (13131) and (I314p describe the same physics they 
suggest (ignoring the signs and numerical factor) the following qualitative 
identifications 

J ~ pI^ (315) 

The second of this equation may seem strange since A (and the covariant 
derivatives) are vectors. The combination that seems to be plausible is 

Di(f)j = di(j)j - eikip(j)k 

There is some confirmation from the Abelian version in the Maxwell-Higgs 
case at self-duality [IQ]. Eqs. fl315p further suggest that the magnetic field 
of the model can be associated with the physical velocity, B ~ v (however 
in this framework no connection can be made with the Elsasser variables 
u = V + B, w = V — B). It is interesting to remark that the magnetic field 
B can as well be associated with the physical vorticity, since the Chern- 
Simons Lagrangean has the unique property that connects directly the field 
tensor F^^, with the current J'^, as is shown by the second equation of mo- 
tion, —ne^^PPyp = iJ^. Then the relationship which is fundamental for the 
connection between the field-theoretical framework and the physical model, 
Inp = ifj (assumed previously as a simple change of variables) appears now 
consistent with the physical meaning of B, since 

5~Alnp = A^ = tu (316) 

This also suggests 

kB ~ pIu (317) 

The detailed form of these identification cannot be made more precise and we 
limit ourselves to a dimensional analysis. In these relationships there is no 
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factor of dimensionality to intermediate between the two sides. The dimen- 
sional factor that multiply the first relation in Eq. fl315p must also multiply 
the second relation, due to Eq.f l314p . As will be verified below, the factor 
(we note it x) must have the dimension 

[X] = L' (318) 
and this implies for the dimansions of the variables {L is length, T is time) 

[X][J] = [pI^]=y (319) 
L'[J] = Y 

or 

b1 = ^ (320) 
In the second relation of Eq. (13151) we have 

[x][v'A] = 0;] = ^ (321) 
L'FA] = ^ 

As we have mentioned before, the quantity v"^ is related with the physical 
background of vorticity generated by the gyration of the particles. Then its 
dimension is 

[v'] = ^ (322) 

F] [A] = ^ (323) 
[A] = 1 (324) 



from which we derive 



or 



This further gives 



L2 

All dimensions become coherent if we identify 



[B] = ^ (325) 



n = Cs (326) 
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For example, using again the unique dimensional coefficient Xi Eg- (13171) is 
dimensionally correct. 



m^m = [pIu] (327) 

TL2 T 

One can now verify that all equations in the field model have coherent di- 
mensions. 

We have now a qualitative association between the physical variables and 
the field model variables and we also have the physical dimensions of the 
latter. We note that the covariant derivatives (having dimension L~^) 

D^ = d^ + [A^,] (328) 

cannot have a clear identification in terms of physical variables. One can 
only say that the zero component is 

where we have taken into account the second relation from Eq.f l315p and 
included the unknown dimensional factor L^. From the Eqs.f lll2"i ) and f ll3ip . 
(11321) we note that it is not possible to express in terms of the classical {ipy 
V, u) variables the potentials and Ay. 



10.5 Comment on the physical constants and normal- 
isations 

One of the characteristics of the physical model is the presence of a uni- 
form background of vorticity. In the absence of any excitation we have on 
any contour in plane a tangential projection of the velocity of the particles 
performing the Larmor gyration. 

An arbitrary contour (say, a large circle of radius R) will intersect the 
circle of the Larmor gyration (of radius ps) and one can calculate an average of 
projection of the velocity onto the tangent at the contour line. Supposing that 
-R ^ Ps, the contour intercepted by the Larmor circle can be approximated 
with a stright line that intersects the circle between the angles and tt — 6q. 
The contour is a chord and the average of the velocity's projection on it. 
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Vc {0), is 



° d9 

""o = / (330) 



1 



7T-29o 

71-200 



dOvL sin 9 

00 



COS ^0 



where vl is the velocity on the Larmor circle. Now we can average over the 
various lengths of the chord inside the Larmor circle, 

V ^ ^ COS 6^0 (331) 

Jo TT TT - 20o 

2vl r'"^ sinr , 
ar 



TT 

2 (1.37) 



TT 



-vl 



The symmetric situation will bring a similar factor and finally the average 
projected velocity is within a factor not far from unity equal to vl- Now 
consider the definition of the rotational 

a; = V X V = lim^^^- (332 

where A is the area inside the closed contour F. We have, within a unity-size 
factor 

V ■ dl ^2tiRvl (333) 
A = nR^ (334) 

u ~ lim — (335) 

Vl = Ps^ (336) 



Then 

Since 
we obtain 



a; ~ Q lim ^ (337) 
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i.e. we obtain that the value of the vorticity in a region with uniform density 
of Larmor gyrating particles is Q, the cyclotronic velocity. 

We have in this moment three parallel models, representing the same real- 
ity, which we call the Charney-Hasegawa-Mima vortical flow. The connection 
between these three models implies a comparison of the physical quantities 
present in each of them. For this reason we have to consider the physical 
content of the field-theoretical model and in particular we will introduce 
nondimensional variables. The two physical quantities appearing explicitly 
in the field-theoretical model are k, and f^. 

10.6 Comparison with numerical simulation and with 
experiment 

The second factor of the nonlinearity i.e. (cosh?/; — p) in all versions of the 
equation derived above, in particular in Eq. (12171) can also be negative, under 
a certain choice of normalizations. Then, a certain aspect of the graphs 
resulting from numerical simulations (see Seyler [lO]), i.e. the presence of 
two visible symmetric extrema on the graph {vorticity, streamfunction) is in 
agreement with our equation. 




Figure 6: The nonlinear term in the equation, for p = 7 and p = 10 



It is interesting to remark that similar pictures to our result have been 
found in the study of geostrophic turbulence [H]. The scatter plots for the 
pair {oOjip) obtained from experimental study of decaying vorticity field and 
represented in the figure 21 of this reference are very similar to our result 
(with the choice of inverse sign for the streamfunction). 
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10.7 The vacuum fields 

We comment on the meaning of the vacuum value of the matter field. Ob- 
viously, it must be related with the presence, in the physical model, of a 
background of vorticity simply given by the gyration of the particles. This 
is the sense of the fact that, even at infinity, we have a constant density of 
matter, |0|^ = f^. The excitations in the form of large vortices take place 
on this background, whose value of vorticity is very high, the ion cyclotronic 
frequency, flci- As we said before we can possibly identify 

= rici (338) 

K = Cg 

The physical vorticity is derived from 

F+- = 1 - 2 (pi + p,)) (pi - p2) H (339) 

10.8 The subset of self-dual states of the physical sys- 
tem 

We should remind that the identification 

^Z-^lnp (340) 

was done after the equations of motion in the field theoretical framework 
have been reduced to the equations of self-duality and stationarity. Therefore 
it is not surprising that, returning with Eq. fl340p to the original, Charney- 
Hasegawa-Mima equation, we find that this one is verified by the functions 
p obeying our equation fll44p or Eq. 01481) for 

[(-W X fi) • V] V^V' = (341) 

because 

[{—WiIj X n) ■ V] [— sinh il) (cosh — 1)] = 

The fact that, for any solution if) of the Eq. (11481) the equation of Charney- 
Hasegawa-Mima is verified at stationarity is useful as a confirmation but is 
of moderate significance, due to the large space of functions that can verify 
Eq. fl34ip . The subset of self-dual states is much smaller and precisely defined 
by Eq.dUHD. 
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10.9 Comment on the self-duality 

In general the self-duality should be seen as a property of a particular geometrico- 
algebraic object : a fiber bundle. This consists of a basis manifold (on which 
one has to define an atlas of compatible charts), a fiber attached to every 
point of the basis manifold (in physics the fiber is seldom said space of inter- 
nal symmetry) and a group of automorphism of the typical fiber. The local 
structure on the basis and in the fibre space is Euclidean since they both 
are manifolds. One can construct the total space of the fiber bundle, which 
is locally a Cartesian product of an open set of the basis with the space of 
the fiber, and a projection operator acting in this total space and projecting 
the points of the total space onto the basis. The transition functions be- 
tween neighboring charts consist of elements of the group. A connection is 
a differential one-form defined in every point of the total space and taking 
values in the algebra of the group. The curvature is the differential two-form 
obtained by an exterior differentiation of the connection one-form. For a 
concrete example, the connection is the potential and the curvature is 
the field strength F^^ like in the electromagnetism, or in any other theory 
expressed in similar terms. There is a Hodge duality operator, denoted * : 
applied on a differential p-form in a space with n dimensions, it generates 
a differential n — p form, such as the exterior product of these two forms 
produces a scalar multiplying the unique n form that can de defined on the 
n-dimensional space, i.e. a multiple of the volume form. 

The self-duality is the property that consists of the equality between the a 
differential form and its Hodge dual; this naturally requires that the space be 
of even dimension. Only the differential forms of the order representing half 
of the (even) dimension of the space can be self-dual since only in this case 
their duals will be of the same order. For example in a space with dimension 
four, differential two-forms can be self-dual. 

In particular physical models the differential two-form representing the 
field strength and its dual are equal at self-duality. But this two-form repre- 
sents the curvature of the fibre bundle, therefore at self-duality the curvature 
is equal to its dual. In many cases, this equality is realised by the fact that 
the curvature is zero and one says that the space is flat. 

When the self-duality is realised as a condition of flatness it is possible 
to express this equality as the compatibility condition of a system of linear 
differential equations. This makes possible to introduce a Lax operator and 
the self-duality equation is exactly integrable by Inverse Scattering Trans- 
form. A set of infinite invariants can be found. One example of this type is 
the classical sigma model. 

In particular cases, (hke ours) the geometrical structure is less clear. The 



74 



self-duality is expressed by the fact that the action functional is minimized 
i.e. the Bogomolnyi limit is saturated. 

10.10 Comment on the 6*^ order potential 

The Abelian version of this theory but with the Maxwell term instead of 
Chern-Simons is well known from superconductivity theory. It implies a 
potential of only fourth order which provides the symmetrical vacua of the 
theory and allows mass generation for the Maxwell photon via the Higgs 
mechanism. 

However in the present theory a sixth order potential is necessary. It 
has been demonstrated that with only a sixth order potential one can have 
self-dual states. This has been shown by a simple verification which we 
have reproduced in the Section about the energy functional related to the 
Lagrangean density. 

However the necessity to include a sixth order potential in the Lagrangean 
density has a profond origin. This has been shown in series of papers [5U] . 



It has been shown that the Bogomolnyi lower bound for the energy and 
the first-order-in-time differential equations obtained at self-duality are a 
property of a classical field theory which possesses a topological charge. The 
theory is a reduction of a supersymmetric (susy) theory in which the topo- 
logical charge appears as the central charge of the susy algebraic structure. 
( A susy theory is a classical field theory in which besides the usual fields 
there are other field- variables with the property that they anti-commute, i.e. 
they are classical spinors). It is interesting the way in which this has been 
shown [SI],[S2]- First it is shown that any supersymmetric theory which 
possesses a topological charge necessarly possesses a Bogomolnyi bound and 
SD equations of motion. Then for a given field theory where a topological 
conservation charge exists, it is first constructed a supersymmetric extension, 
adding the anti-commuting variables and other variables that are necessary 
to close the new algebraic structure. In this extended theory the central 
charge of the susy algebra is the topological charge of the initial theory. The 
Bogomolnyi bound is identified. Finally it is shown that returning back from 
the susy extension to the original theory, one still preserves the Bogomolnyi 
bound. The relation between the potential W in the extended theory and 
the potential U in the classical non-susy theory is 



m 




(342) 
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and for a symmetrical two- vacua potential W we have a potential U of sixth 
degree in 

Lee, Lee and Weinberg [53] show explicitly in an Abelian case how this 
form of the potential is obtained from the requirement that the model can be 
extended to a = 2 supersymmetric model. They begin by constructing an 
= 1 supersymmetric generalization of this Chern-Simons Higgs theory, in 
which the form of the potential f [\<f)f) is not yet specified. Adding a single 
pair of Grassmannian varibales (6', 6') to the set of variables of the original 
model requires to extend the model by introducing additional fields. This is 
necessary since the supersymmetry transformation must be closed and the 
original fields are not sufficient. The model will include a matter super-field 
$ which consists of: a complex scalar field 0, a complex spinor field ijj and 
an auxiliary scalar field F; a real spinor field which contains a real photon 
field Afj_ and a Majorana spinor photino field A. The action is 

5 = y d'^xdt j dedel^-^KD'^T'^DpVa (343) 

+/($*$)} 

where the first term is the generalization of the Chern-Simons term. The 
integration over the Grassmann variables 9 and 9 can be done explicitely 
and in this process it is required to make them visible in the expression of 
the superfield In this way there appears in the action density first and 
second order derivatives of the potential function / since only in this way 
(by this Taylor expansion) the Grassmann variables will appear explicitely 
and can be integrated. The action becomes 

S = j d'xdt!^^Ke''''^A^F,x + iD''<P)iD^(P) (344) 
-^kAA + i^YDf,ip + i (^A0 - Xip(j)*) + F*F 
+f {F*<p + F<p*) - \f" (<p'i^r + 

where the superscript c means charge conjugate and the prime means deriva- 
tive of the function / to its argument, The equations of motion for the 
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fields A, F and F* are 



K 



(345) 



A 

F = 
F* = -07' 

They permit the replacement of the corresponding functions in the action. 



S = j £xdt!^^K€^''^A^F,x + iD^'(P)iD^(P) 



(346) 



+ 



,2 / 1 



--/" -/' 



In order this action to be invariant under an = 2 extended supersymmetry 
it is required that the term on the third line in the above formula vanishes 

/" = -- (347) 



K 



or, integrating two times on the variable ^ = |$| 



2\2 



(348) 



1 

'2k 



where v"^ is a constant. The form of the potential term in the action is then 
obtained from the first term in the second line of Eg. (13461) 



1 2 r,2 



2k, 



-V ) 



(349) 



The action contains a bosonic part that has this potential and this is actually 
the Abelian version of the model discussed in this work. 

The fact that at self-duality the theory is a part of a larger super symmetric 
theory and that this explains the form of the scalar self-interaction may help 
us to trace the meaning of the changes we find between the sm/i-Poisson 
equation (for the ideal fiuid) and the double- sinh-Poisson equation, for the 
fluid of ions with Larmor gyration. 

The special Higgs potential that appears in the Lagrangean density and 
leads to self-dual states has a symmetric minimum which is degenerate with 
the symmetry-breaking one. This means that the system can have non- 
topological solitons which verify the same self-dual equations. 
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11 Appendix A : Derivation of the equation 

Consider the equations for the ITG model in two-dimensions with adiabatic 
electrons: 

^ + V-(v,n,) = (A.l) 

^ + (v,-V)v, = A(_v0) + Av,xB 
at rrii rrii 

We assume the quasineutrality 

Hi He (A. 2) 

and the Boltzmann distribution of the electrons along the magnetic field line 

He = no exp f -li^ ) (A.3) 



The velocity of the ion fluid is perpendicular on the magnetic field and is 
composed of the diamagnetic, electric and polarization drift terms 

Vi = (A.4) 
T, 1 dn,^ 



\e\ B Hi dr ^ 

— Vo X n 



+ 



B 

1 f d 



B^i \dt 

The diamagnetic velocity will be neglected. Introducing this velocity into the 
continuity equation, one obtains an equation for the electrostatic potential 
<f>. 

Before writting this equation we introduce new dimensional units for the 
variables. 

\p\ fkphys 

^hys (A.5) 

(rpPhys „.phys \ 
, (A.6) 
Ps Ps J 

tphys ^t' ^ t^^y'^i (A.7) 
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The new variables (t, x, y) and the function cf) are non-dimensional. In the 
following the primes are not written. With these variables the equation 
obtained is 

I (1 - Vi) (A.8) 

— ( — V_L0 X n) ■ KBr 

- [(-Vx0 X fi) ■ V J Vl0 
= 

where 

/te,. = — V_L Inrio (A. 9) 

(|57j). Before continuing we compare this equation with the equation of 
paper [58], Eq.(16). Here taking still the units to be physical, the form of 
the latter equation is (Eq.(12) from that paper) 



d \e\(j) d 1 



dt Te dtBVti 

H ^ V_L In no 

— V_L</) X n ^ |e| 



(A.IO) 



B " Te 

1 







The term containing the gradient of the equilibrium density comes from the 
continuity equation, as convection of the equilibrium density by the fluctu- 
ating E X B velocity. The adiabaticity has been assumed. 



n 



no Te 



{A.ll] 



and we consider that the temperature is constant (the calculations can easily 
include a dependence Te (x) ) . 
For the second term we have 

^ ^9 , 1 Tg lei (/) 11 To lei 

Vi0 = i^^Vl^ = --^-^Vl^ (A.12) 



BQi BVLi\e\ ^ T^ a, ^ m,; ^T 

1 



rrii 
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This will become (with its sign) 



dt 



7'2 ll 



(A.13) 



in the new variables Eqs. flA.5p - flA.7p . 
The third term is 



X n 



B 



V_L In no 



B \e\ 
1 



X n ) ■ V± Inno (A.14) 



\e\B 
mi 

VL — 



X ii I ■ V_L In no 



T 

|e| (/) ^ , 
-Vx^^ X n ■ Vxlnno 



X n I ■ (-V^lnno) 



This will become 



- Qi (-ps V_L0 X n) ■ (-p^ V_L In no) (A. 15) 

and in the normalised space variables 

- a (-V10' X n) ■ (-Vllnno) (A.16) 
The last term (with the polarization nonlinearity) is in physical units 

- [(-Vx0 X n) . V J Vl0 
This is converted to non-dimensional variables 



1 



Ps |e| 



fPsVj 



T, 



X n 



ps 



-PsVj 



Collecting the physical coefficient we have 



T 



111 

pI |e| 



1 fT, 



-P^vl 



|e|0 



Pf 



Pt 



(A.17) 



(A.18) 



(A.19) 
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Then, in the normahsed variables, this term becomes 

X n) ■ V'^] {-Vl<P') (A.20) 

Then the Eqs. f[O0|) with the new form of its terms (lA.lSp . (1A.16P and 
(lA.20p becomes 

-n, (-Vi0' X n) ■ (-V^lnno) 
+n, [(-Vl</)' X n) . VI] (-V'i0') 
= 

Introducing the time unit and ehminating the primes 

I (1 - Vi) (A.22) 

- (-V±0 X n) ■ (-Vxlnno) 

- [(-Vx0 X n) ■ VJ Vi0 
= 

The last term is the convection of the vorticity 

cu = Vi0n (A.23) 

by the velocity field 

\E = -V_L0 X n (A.24) 
We use the definition Ke^ = V _l In no or 

Key = — n X Vxlnno (A. 25) 

Then the resulting equation is 

(1 -Vl)^-K^- [(-Vx0 X n) . VJ Vi0 = (A.26) 

The same equation but without the linear (density gradient) term has 
been derived as the "shielded convective ion cells" [TT] 

(i-pM)f -[0,pM0]=o 

and as a possibility to describe the Kelvin-Helmholtz instability modified by 
the finite parallel electric field and its associated current density with 



e- 



dt Te dt 
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where the potential is not normahsed yet. The enstropy is conserved 
[/ = J d^r [{psV±(l)f + (p'Ax0)'] = const 

12 Appendix B : The Euler-Lagrange equa- 
tions 

The calculations from this Appendix should be considered as a guide for 
a first contact with the methods of the theory of non-Abclian gauge field 
interacting with nonlinear (= self-interacting) scalar matter field. Here the 
calculations arc not pedagogical (in particular we treat asymmetrically the 
fields Afj^ and AjJ and we suggest that after these first steps other lectures 
are necessary, from field-theory genuine sources. 

Consider again the Lagrangean density 

C = -ne^''Ptr(^d^A,A,+ '^A^A,Ap^ (B.l) 

-tr [(D'^<l>)^ (L>^0)] 
-V (0, </.t) 

The functional variables are 

Ao,AlA,,AlA,,Al (B.2) 

and they are all SU (2) matrices with complex entries. 

12.1 The contributions to the Lagrangean 

12.1.1 The Chern-Simons term as a differential three-form and 

the presence of a metric 

Apart from a factor, the gauge Lagrangean is the trace of the Chern-Simons 
differential three- form on a principal bundle with group SU (2). 

n = -^tr ( AAdA-^A A AA A ) (B.3) 
oTT^ y 3 J 

1 ^^.up ( A^d^A^ - -e^^^AlAlA^] d^x 



167r2 
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The trace of the Chern-Simons form can also be expressed using the 
exterior differentation and exterior product of forms 



cs{A) = — [ tr ( A A dA + -A A A A A ) (B.4) 
4vr V 3 y 

where, for three algebra-valued differential one-forms Ak, k = 1, ...,3, one 
has 

tr {A, AA2A As) =^ ^tr {A, A [A2, A3]) = ^tr {[A,, A2] A A3) (B.5) 

A factor of 1/2 can also be extracted from the first term in flB.4p if we add 
minus its expression but with two of the three indices exchanged. Then 



cs 



;A) =i- J E'^-'ftT (^A^ {d^A^ - dpA,) + [A,, A;^ (B.6) 



The normalizing constant in Eq. (lB.6p is related with the fact that the 
integral of the Chern-Simons form is a topological invariant for adequate 
boundary conditions and has integer values. For what we need, the gauge 
field Lagrangean can be taken such as to lead to the gauge-part in the action 



S, = JdtJ £x j-i/t^^-^tr (^A^ id,A, - d,A,) + [A,, A^]^ | (B.7) 
Therefore we will use the following expression 

= —ne^'^'ti (^A^ {d^Ap - d,A,) + '^A^ [A,, A^]^ (B.8) 

We write in detail Eq.f lB.Sj) . The first term is 

e'^-'P [A, id,A, - d,A,)] (B.9) 
= e'''Ao{d,A2-d2A,) 
+e'''Ao id2A,-d,A2) 
+e'^^A, idoA2 - d2Ao) 
+e'^'A, id2Ao - doA2) 
+e'''A2 id,A, - doA,) 
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or 

e^''P[A^{d,A,-dpA,)] (B.IO) 

-Ao (92^1 - diA^) 

-Ai (^0^2 - ^2^0) 
+Ai (^2^0 - ^0^2) 

-^2 {diAo - d^Ai) 
+^2 {doA^ - aiAo) 

This is simply two times every distinct term in the sum 

e^^^p[A^{d,Ap-dpA,)] (B.ll) 
= 2A^{diA2)-2A^{d2Ai)-2Ai{d^A2) 
+2Ai (^2^0) - 2A2 {diAo) + 2A2 {doAi) 

We continue by calculating the second term in the CS action 

e^^PA^[A,,A,] = e^^^PA^{A,Ap-A,A,) (B.12) 

(A0A2 - A2A0) 

+6^20^^ (^^^^ _ 
+£210^2 (AiAo - A,A,) 



or 

= (B.13) 

= Aq{AiA2- A2A1) 

-Ao {A2A1 - A1A2) 
-^1 {A0A2 - A2AQ) 
+Ai {A2AQ - AqA2) 
-A2 {A^Aq - AqA^) 
+A2 {AqA^ - A^Aq) 

This is actually two times every distinct term in the sum 

]^e^^-''A^[A,,A,] = A^A,A2-A,A2A, (B.14) 

-A1A0A2 + A1A2A0 
-A2A1A0 + A2A0A1 
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It results that the form of definition 



£i = -^Ke^^^^tr (^A^ (d^Ap - d,A,) + ^A^ [A,, A,]j (B.15) 
can be written 

£i = -KtT{Ao{d,A2)-Ao{d2A,)-A,{doA2) (B.16) 

+A, {d2Ao) - A2 (diAo) + A2 (doA,) 

2 2 2 

+-A0A1A2 - -A0A2A1 - -A1A0A2 

2 2 2 

+ -A1A2A0 - -A2A1A0 + -A2A0A1 

We recall that every function An is actually a matrix, in the adjoint repre- 
sentation of the SU (2) algebra. We can use the property of invariance of 
the operator Trace of a product of matrices to a cyclic permutation of the 
factors of this product. The third term in the first line without derivatives 

A1A0A2 A2A1A0 A0A2A1 (B.17) 

The second term in the second line without derivative 

A2AiAo A0A2A1 (B.18) 

These two terms will add to the second term of the first line, giving 

- '^AoA2A, - ^-A^A^A2 - ^AsAiAo ^ -2A0A2A1 (B.19) 

The first term on the second line without derivatives 

^lAsAo ^ (B.20) 

The last term 

^ ^ A0A1A2 (B.21) 

These two terms are added to the first term of the first fine without deriva- 
tives 

2 2 2 

-A0A1A2 + -AiA2^o + -^^0^1 ^ 2^o^i^2 (B.22) 
Finnaly, we collect all terms that do not contain derivatives 

^^^1^2 - ^AoA2Ai - ^A^^ (B.23) 

2 2 2 

- -A2A1A0 + -A2AqA^ 

= 2AqAxA2-2AqA2Ax 
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At this point the gauge-field Lagrangean is 



A = -KtT{A,{d^A2)-A^{d2A,)-A,{d,A^) (B.24) 

+Ai (^2^0) - A2 (9iAo) + A2 (9oAi) 

+2A0A1A2 - 2AQA2A1} 

Change of the form of the gauge-field part of the Lagrangean, from 
integration by parts Since the Lagrangian density is integrated in order 
to obtain the action functional, we can consider the effect of integrations 
by parts. These operations will move the differential operators between the 
factors of the monomials appearing in the expression of the Lagrangian den- 
sity and will also generate boundary terms. In general the boundary terms 
are zero for well-behaved functions, but in our case the presence of a fi- 
nite condensate of vorticity at infinity can produce finite terms. We will 
develop below a calculation based on integration by parts, removing the spa- 
tial derivatives from acting upon Aq but we will ignore the boundary finite 
terms. Then the calculation is simply useful for the comparison with other, 
well known, forms of the CS Lagrangian. We will not use the form of the 
Lagrangian derived from these operations Eq. flB.28p . for obtaining the Euler- 



Lagrange equations, and we will rely on Eq. flB.24p . 



We turn to the terms containing derivatives. It is possible to make a 
integrations by parts using the formula 

d , , ^dZ dY* 

— YZ =Y*— + -— Z B.25 

dx dx dx 

We apply this for the following two terms and furthermore we use the cyclic 
symmetry inside the Trace operator 

A,{d2Ao) ^ d2{AlAo)-{d2A,)Ao (B.26) 
^ -Aoid2A,) 



-A2{diAo) -di{A;Ao) + {diA2)Ao (B.27) 
^ AoidiA2) 

We collect all terms containing derivatives 

Ao {d,A2) - Ao {d2Ai) - Ai {doA2) 
+Ai id2Ao) - A2 id,Ao) + A2 (9oAi) 
= 2Ao idiA2) - 2Ao (Sa^i) - A^ (^0^2) + A2 (doA^) 
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Finally, the gauge-field Lagrangean density results 

A = -KtT{2Ao{diA2)-2Ao{d2A,) (B.28) 
-A^ idoA,) + A2 (doA^) 
+2AoAiA2 - 2A0A2A1} 

Comparison with known forms of the gauge Lagrangean Let us 

check this formula by comparing with the situations where the space part is 
separated [i9] . 

The metric is defined from the general expression of the differential length 
in the 2 + 1 dimensional space, as 

c/s2 ^ _^^2 ^ h..dx'dx^ (B.29) 

{h}^ is the space-part of g^") and one can separate the spatial and temporal 
parts of the action. In our case the metric is g = diag {—1, 1, 1}. The action 
is 

5i = j dt j d'^ xti {~Ke'^AidoAj + ne'^AoFij] (B.30) 



where 



and 



= d,A, - djAi + [A, Aj] (B.31) 
F,o = D,Ao - doA, (B.32) 



Then, in the integrand 



^21 

From Eq. (]B.33p where h is the metric, we have 



e'^ = (B.33) 



-e'^AidoA, + e'^AoFij (B.34) 
-e'^A,doA2 - e^'A2doA, 
+e'^AoFu + e^^AoF2 



Now we replace 



-e'^AidoAj + e'^AoFij (B.35) 
-AidoA2 + A2doAi 
+A0F12 — A0F21 



Fij = diA,-djAi + [A,Aj] (B.36) 
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and obtain 



+Ao (91^2 - 92^1 + A1A2 - A2A1) 

-Aq {d2Ai - diA2 + A2A1 - A1A2) 



(B.37) 



We note that some terms are repeted 



-e'^AidoA^ + e'^A^Fij 

-AidQA2 + As^oAi + 2AQdiA2 - 2Aod2A 

+2A0A1A2 - 2A0A2A1 



(B.38) 



The result is identical to Eq. (]B.28p 

£1 = -Kti{-A^doA2 + A2doAi + 2AodiA2-2Aod2A^ (B.39) 
+2A0A1A2 - 2A0A2A1} 

We also note that until now there was no need to consider summation 
over components of vectors using the metric coefficients. In case of a product 
of the form x^x^ it will have to consider the metric. 

12.2 The matter Lagrangean 

The form is 



C2 = -tr [{Dyy 
Using Eqs. flH^ . flHTl) and flHUj) we can calculate in detail 



(B.40) 



£2 = -tr {D^'4>y {D,4>) 



(B.41) 



= -tr 
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We have to expand the products 



at 



dx dx dx dx 
ox 
ax 



9y 9y 9?/ dy 
dy 
9y 

12.3 The Euler-Lagrange equations 

The Euler-Lagrange equations 

9 SC SC 



(B.43) 



We use distinct notations for the three components of the Lagrangean density, 
C = Ci + C2 — V where £1 is the gauge field part, C2 is the "matter" part 
and V is the nonlinear self-interaction potential for the "matter" field. We 
use the detailed expressions for Ci from Eq. (]B.24ll and £2 is given by the 
Eq. (IB.42p . The functional derivations are done separately on these two parts. 



12.3.1 The formulas for derivation of the Trace of a product of 
matrices 



Use the formulas (see Ref. [i7] ' 



Atr (AX) = (B.44) 
Atr (XA) = A^ (B.45) 
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-^tr (X^A) = A (B.46) 

Atr (AX^) = A (B.47) 

-^tr (AXB) = A^B'^ (B.48) 
uX 

— tr (BX^A) = AB 
rfX ^ ^ 

-^tr (XAX^) = X (A + A^) (B.49) 
ctX 

-^tr (X^AX) = (A + A^) X (B.50) 
dX 

-^tr (AXBX) = A^X^B^ + B^X^A^ (B.51) 
uX 

-^tr (AXBX^C) = A'^C^XB^ + CAXB (B.52) 
ctX 

where A, B, C, X are arbitrary complex matrices. 

12.4 The Euler-Lagrange equations for the gauge field 
12.4.1 The variation to Aq 

The equation of motion resuhing from the variation to Aq is 

' (B.53) 



dx^6{^) 6Ao 



dx^ 6 (doAo) ^ dx^ 6 {diAo) ^ dx^ 6 (SzAq) 'sTo ~ ^ ' ' 

Functional derivatives at Aq of the gauge field Lagrangean The 

gauge field Lagrangean is Eq. flB.24p 



A = -KtT{Ao{d,A2)-Ao{d2A^)-A^{doA2) (B.55) 
+A, (d^Ao) - A^ id,Ao) + A^ (doA,) 
+2A0A1A2 - 2A0A2A1} 



and we have to calculate 



dx0 6{doAo) dx^5{diAo) dx^ 6 {d2Ao) 6A0 
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The first term in the Euler-Lagrange equation flB.56p for Aq is zero since 
For the second term there is only one contribution 

= (-K)(-a,Aj) 

The third term also consists of one contribution 
d 6Ci 8 6 



-K)ti{A,{d2Ao)} (B.58) 



= {-n){d2A^) 

The last term in Eq.(??) is the derivative of Ci to the functional variable 

^ = -K^tT{Ao{^^A2)-Aoi^2A^)-A,{^oA2) (B.59) 
dAo dAo 

+A^ id2Ao) - A2 i^^Ao) + A2 (doA,) 
+2A0A1A2 - 2A0A2A1} 



In detail, every term 



^ tT{Ao{d,A2)} = id,A2f (B.60) 



SAo 

tI{-Aoi^2A^)} = -i^2A,Y (B.61) 

tT{-A,{doA2)} = (B.62) 

tr{Ai(92Ao)} = (B.63) 



5^0 

6 



SAo 
6 



SAo 
6 



tT{-A2{diAo)} = (B.64) 



5^0 

91 



Mo 



ii{A2 (9oAi)} = 



5A, 



^ ti{2AoA^A2} = 2{A^A2Y 



6 



SAo 



ti{-2AoA2Ai} = -2 iA2A^] 



(B.65) 
(B.66) 
(B.67) 



Collecting these formulas we find 
6Ci 



5Ao 



i-K) {{d^A^f - (d^A.f + 2 {A,A2f - 2 (^2^1)^} (B.68) 

{-k) {diA2 - diAi + A1A2 - A2A1 
A1A2 - ^2^1}^ 

(-/€) (9iA2-92Ai+[Ai,A2]f 
+ (-K)([Ai,A2]f 



The result of the variation of the gauge-part of the Lagrangian £1, to the 
functional variable Aq is obtained from the results Eqs. (]B.57l) . (1B.58I) and 



d 5Ci d 6Ci 

+ 



5Ci 



(B.69) 



dx^S{diAo) dx^6{d2Ao) 6A0 

~{-k) {d,A2~d2A, + [AuA2]f 
-i-K)i[A^,A2]f 

= Kid,A2-d2A, + [A,,A2]f 

+ i-K) [-diA^ + d2Al - [Ai, ^2]^} 

= 2K{diA2-d2A^ + [A^,A2\f 

= 2KiF,2f 

Functional derivative with respect to Aq of the "matter" Lagrangean 

We continue with the variation to the functional variable Aq of the of the 
"matter" part of the Lagrangean £2 is 



d 6C2 d 6C2 d 6C2 



6C, 



dx^didoAo) dx^6{diAo) dx^ 6 {d2Ao) 6A0 



(B.70) 



where 



£2 



-tr 



92 



has the detailed expression given in Eq. flB.42p . The Lagrangean is 
£2 = -tr 



dd) dd)^ , , dd)"^ , , 



dt dt dt 



dt 



dt 

d(f)^ dcj) d(f)^ 



+ 



90t 



dx dx dx 



dx 



Ml 



ox 



+ + ^^20 - ^0^2 



5|/ dy dy 



dy 



+ 



dy 



5y J 



Calculation of the variation of the matter Lagrangian to the field 

Aq The first term is 

since there is no explicit dependence of £2 on (doAo). 
The next two terms in the variation of C2 are 

d 6C2 



dx^ S (diAo) 
d 6C2 
dx^ S {d2Ao) 








Again, there is no dependence of £2 with respect to diAo and c?2Ao and these 
contributions are zero. 
The last term is 



6C2 

SAo 



6An 



-tr 



iD^<l>y {D,<P) 



(B.71) 
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Only few terms from have non-zero contributions 
We calculate in detail every term 



5Ao 1 dt ''^ \ \ dt 



T 



tTl-^Ao<p[= i^] {<pf (B.73) 



^ -tri^Mo\ = -{^<^] (B.74) 



5Ao [ dt^ \dt 

- ^tr {0tA°tAo0} = - (0tA°t)^ (0)^ (B.75) 

In this is formula we have applied Eq. flB.48p with A A^'^ , B =0 since the 
functional variables ^4°^ and Aq are independent. 

- T^tr {-0tAOt0Ao} = (B.76) 

OAq 

In this is formula we have applied Eg. (IB. 451) with A =(f)^A'^'^(l), as explained 
above. 

- ^tr {-A^^<p^Ao<P} = (AOt0t)^ (0)^ (B.77) 
This is formula flRlHl) with A =A°t0t^ b =0. 

5 



5/1 



-tr = - (B.78) 





Here the Eq. (IB.45l) has been used with A = Aq(^ 

Now we sum the results from Eqs.( lB.7^ to (IB. 781) 



6An \ dt J \ dt 



T 
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From the first, third and fifth terms we separate to the right the factor (0)^, 
and similarly in the other terms. 



6C2 ^ \fd^ 
SAo \\dt 



'dt 



Now, taking the transpose operator out of the paranthesis. 



5C2 



dt 

a0t 



(B.81) 



5C2 
Mo 



~dt 



(B.82) 



+ 







We now change the upper index into the low index for the potential A\ 



and obtain 
5C2 



5A^ 



+ 



84)^ 



dt 



+ - At0t 



(B.83) 



at 







+ 



- (Ao</>)^ + {<j>A,)^ 



dt 



dt 

■ + ^00 - (f>Ao 







•_ 









at 



+ Ao(l) - Mo 



T 



(D 



{0(L>o0)^}^-{(L>o0)V} 
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Collecting these results 



d 5C2 d 5C2 d 5C2 



5C2 



dx^d^dQAo) dx^d^diAo) dx^5{d2AQ) 5Aq 



;b.84) 



0, Po0)^]} 



T 



Calculation of the variation with respect to Aq of the scalar poten- 
tial This is extremely simple since the scalar potential does not depend on 
Aq nor of its derivatives to x^. 



5V 







(B.85) 



Final result for the variation of the Lagrangian with respect to Aq 

We assemble the partial results: Eq. (]B.69p and Eq. (lB.84p 



8 5 8 5 8 5 



8x^5{8qAo) 8x^5{8iAq) 8x^5{82Aq) 5Aq 



(£1 + C2 (Blg^) 



The equation is 



2K{F^2y 



or 



0, PO0)^ 



2kFi2 

The left hand side can be written 

and we change the order of the terms in the commutator 



KB 



up 



po0)^0 



(B.87) 
(B.88) 

(B.89) 
(B.90) 



and multiplying with — 1, 



Oup p 



up 



{B.91] 
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We note that in the right hand side we have a part of the expression of the 
current 

Jo~-^{-[W)^0]} (B.92) 
according to the definition of the current 

-Ke^'^pp^p = -iJo (B.93) 

Then at this point of the derivation it is suggested the following form of the 
= component of the equation of motion 

- Ke^'^'PF^p = i J° (B.94) 

Below we will join to this part of functional variation another part, re- 
sulting from the functional varaition with respect to A^^ . 



12.4.2 The functional veiriation with respect to the Vciriable ^4°^ 

We have to calculate 

8 5 8 5 8 5 8 \ ,n . .r^ 



8x^5{8qA^^) 8x^5{8^A^^) 8x^5{82A^^) 5A^^ 

(B.95) 

Each term is calculated separately. 

Functional derivatives of the gauge field Lagrangian with respect 

to ^4°^ The gauge-field part of the Lagrangian is 

A = -«tr{Ao(aiA2)-Ao(92A)-A(9o>l2) (B.96) 
+^1 (92A0) - A2 {d^Ao) + A2 (9oAi) 
+2AqAiA2 - 2AQA2A1} 

We note that the gauge field Lagrangian Ci is not expressed in terms of 

f 8 5 8 5 8 5 



^axO(5(aoAot) dx^5{diA^^) dx'^5{d2A^^) 5A^^ 



97 



Functional derivatives of the "matter" Lagrangian with respect to 

y^of Pqj, ^Yie matter part of the Lagrangian we have to calculate 

8 6 8 6 8 6 6 , ^ 

+ ^ . ..^ .... + ^ . r.r. .... - 1-7^ A 



8x0 6{8oAo^) 8x^6{8iAo^) 8x^ 6 {82A0^) 

where £2 is given in Eq.( ]B.42"ll . 

We have 

-C2 = (B.98) 



8x0 5 [doA^^) 

8 6 

8x^ 6 (Si^ot) 

8 6 



£2 = (B.99) 
£2 = (B.lOO) 



8x^ 6 (aa^ot) 

Few of the terms in Eq. (lB.42p can provide a non-zero contribution 
6C0 



tr|0tA°t^-AVf (B.lOi: 



In detail, the terms are 

6 



6 



-tr 





(B.103) 


(0t)^(Ao0)^ 


(B.104) 




(B.105) 


= {<P^Ao<l>f 


(B.106) 




(B.107) 
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The results are now added 

-(0t)^(Ao# 

In three terms we left-factorize (0^)^ 

M^|-(f)^-(>lo*r+(0>lof} (B.109) 

= w"[-{^y-('i''^'-<i'')] 

= -(,.rll+[A„.^l}' 

For the other three terms the result is similar 

+ (A^^'f - (0*)^ (B.llO) 

= {(f)%i..,.r}(.r 



Then 



5C2 



(B.lli: 



(^00)' , (0^) 

The Euler Lagrange equation is, for A^^ 

d 6C 6C 



(B.112) 



The Euler Lagrange equation in the case of A^'^ is reduced to only the last 
term 

(B.113) 



and is written from Eq. (lB.llip 

SC2 



(B.114) 



12.4.3 The Euler-Lagrange equation derived from functional vari- 
ation to Aq 

We now collect the results of the functional derivatives from both the gauge 
and the matter parts of the Lagrangean, in the Euler-Lagrange equation 
for Aq and add the zero-valued term resulted from the functional variation 
with respect to A^^ . The formulas to be used are Eq. (lB.69l) . Eq. (1B.84l) and 
Eq. flRTIil) 



5Ci SC^ 



12) 



SAo 5Ao 

{D,<py 



SC2 ] 



T 








(B.115) 



or, 



-2«:Fl2 = -[0^A0] - [0,(A0)^ 
and interchanging the factors in the second commutator 



(B.116) 
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The left hand side is the zero component of a tensorial contraction 

- (£°i2Fi2 + £°'^F2i) = - { [(t>\ Do0] - [(^00)^ , 0] } (B.117) 

We will identify the right hand side as the covariant 0-component of a current 
Jo = -2{[0t,Z^o0] - [po0)^0]} (B.118) 

and this equation takes the form of the Gauss law constraint from the main 

text 

- Ks'^P'^F^^ = -iJo = iJ^ (B.119) 



Therefore we conclude that we have derived the component of the equa- 



tion 



12.4.4 The Euler-Lagrange equation from the variation to Ai 

The detailed equation is 

d SC 5C 







B.120) 



and shows a difference compared to the case of Aq: now there is a dependence 
of Ci on the derivatives of the field Ai. 
We have to calculate 



d S 8 5 8 5 



5 



8x^5{8qA{) 8x^5{8iA^) 8x^5{82Ai) 5Ai 



(£1 + £2 - ^) = 



The functional variation with respect to Ai of the gauge-fie\d part 
of Lagrangian This part is 



8 



5 



+ 



8 



5 



+ 



8 



5 



5 



8x^5{8qA{) 8x^5{8iA{) 8x^5{82Ai) 5 A 



£1 (B.121) 



and £1 is given by Eq.( ]B.24j) 



£1 = -Kii{Ao{8iA2)-Ao{82Ai)-Ai{8oA2) 

+Ai (^2^0) - A2 (9iAo) + A2 (9oAi) 

+2A0A1A2 - 2AQA2A1} 
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Term by term 

d 6Ci 



d 



dx^ 6 {doAi 



The second term in Eg. (IB. 121]) 



dx^S{doAi] 
d SLi 



'-K)ti{A2{doA,)} (B.122) 



dx^S{diA^] 
The third term from the Eq. (]B.12ip 

d 6Ci 8 6 



(B.123) 



dx^ 6 (92^1 



dx^ 6 {d2Ai] 

K {d,A^) 



(-K)ti{-Ao{d2A,)} (B.124) 



The last term in Eq. flB.12ip is 

6Ci 6 



6Ai 6Ai 
In detail 



i-K) {-Ai {doA2) + Ai {d2Ao) + 2AoAiA2 - 2A0A2A1} 

(B.125) 

6 



Ml 



Ml 
6 



{-Ai {doA2)} = - {doA2y 



{A^ {d2Ao)} = {d2Aof 



SAi 



{2AoA,A2} = 2AlAl 



5Ai 



{-2A0A2A1} = -2 



It results 



6Ci 



{doA2f + {d2Aof + 2AlAl - 2 (^0^2)^} 



Collecting the results for all the gauge field Lagrangian 



d 



6 



+ 



d 



6 



+ 



d 



6 



dx^6{d(,Ai) dxH{diAi) dx^6{d2Ai) 6A1 

K {doA^) + n {d2Al) - |- (^0^2)'^ + (^2^0)^ + 2AlAl - 2 (A0A2) 



(B.126) 

T 



2k {- [d^Al) + {d2Al) + AlAl - AlAl] 
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We can write 



2k {- {d,Al) + {82 Al) + AlAl - A^A^} 
-2k {doA2 - 62 Ao - A2A0 + AqAs}^ 
-2K{doA2-d2Ao + [Ao,A2]f 
02 



-2kF^ 



The functional variation at Ai of the matter part of the Lagrangean 

The calculations are similar to the previous case for Aq. 
We have to calculate 

d 6 8 6 8 5 6 , ^ 



dx0 5{8oAi) 8x^6{8iAi) 8x^ 6 {82A1) Mi 



where £2 is given in Eq. (lB.42p . 

We have 



8 S , 

C2 = 



8x0 6 {80 Ai J 
8 6 

8x^6{8iAi)^~ 

8 6^ 

-£2 = 



8x^ S {82A 



SC2 s 
6 



6AI 



-tr 



{D^<I))Hd^<P)\ (B.127) 

- (Do0)^ {Do<j)) + {Di<f)) + (^20)^ (1^20) 

6C2 ^ 6_ 

SAi 6A1 



tr|(A0)t(A0)} (B.128) 



-ti {^Ai(j)- ^(f)Ai (B.129) 



Term by term 



6A1 5Ai [ 8x 8x 



6 (84>\ ^ f8cf>^^ ^ 



6A1 \8x^ J \8x^ 



tr<i^A0H-(^) i<Pf (B-130) 
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- ^tr {-^U^V^i} = (B.133) 
- ^tr {-A'^<P^Ai<p} = {A'U^f {<pf (B.134) 

- ^tr {A^VV^i} = - (^^W)^ (B.135) 
Now we collect all these contributions 

As in the case of Aq equation, we right-factorise from the first, third and fifth 
terms 
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In a similar way we have from the second, fourth and sixth terms 

(1^*)^ + {^'A»4>f - {A'UUf (B.138) 

Then, finally 

Adding all contributions 

The total functional variation (gauge and matter parts) 

r d 6 d 6 d 6 5 \ 

\dx^ 6 {doA,) ^ dx^ 6 {diAi) ^ dx^ d {d2Ai) ~ JTj ' 
r d 6 d S d S 6 \ 

^ V^xo 5 (doAi) ^ dx^ 6 (diAi) ^ dx^ 6 (^2^1) ~ STJ ' 



-2kF, 



02 



(B.140) 







Before concluding the calculation for Ai , and having in mind a possible more 
symmetrical form, we study the functional variations to the field ^4^^ 
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12.4.5 Functional variations to the field A^''' 



The equation is 



( aAit N 



= 



(B.141) 



since only from gauge and matter we expect contributions. But we note that 
the gauge part Ci has no dependence on the fields d^A^^ and A^^ 



S{d^A^^) 



= 



(B.142) 



6Ci 



The variation at A\ of the matter pairt of the Lagrangean The 

matter part does not contain any term with dfj^A^^ which means that it also 
cannot contribute to the equation. The only contribution may arise from the 
variation of matter part, C2, to the field 74^''^.The formulas start with 



SC2 



S 



tr 



tr 



(B.143) 



(^00)^ {Do(t>) + (i5l0) + (i?20)^ (i^20)] 



Only the second term must be retained 

SC2 S 



SC2 



Term by term 



tr{(A</>)^ W)} 



(B.144) 



tr U^A^^-^ + - (B.145) 



S 



tr = - 



5 



{4>^fi4>Ai)' 



(B.146) 

(B.147) 
(B.148) 
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6 



S 



tr {-A^^(f>^Ai(f>} = (0^10)^ 



Finally the sum of contributions 



4 ^ -(^'^'(^) 

+(0t)^(0Aif 



The first three terms can be written in the compact form 



- 
- 



T 



= -(0t)^(Di0) 
and the other three terms 



T 



(90 
dx^ 



T 



107 



(B.149) 

(B.150) 
(B.151) 

(B.152) 



(B.154) 



The full result 



5C2 



(B.155) 



I T 



and the contribution of A^'^ to the Euler Lagrange equation has the form 

5C2 



(B.156) 



12.4.6 The final form of the Euler-Lagrange equation derived from 
functional variation to Ai 

We now collect the results of the functional derivatives from both the gauge 
and the matter parts of the Lagrangean, in the Euler-Lagrange equation for 
Ai and add the zero- valued term from the functional derivation to A^^, Eqs. 
(lB.140p . (lB.156p . Therefore we have to combine the following two results 



2kR 



02 



and 



We can subtract the two terms which leads to 



-2kR 



02 



(B.157) 
(B.158) 

(B.159) 



We have 



= 



_102 p I ^120 p 

—Fq2 + F20 
— 2-fo2 



from which 
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We can freely replace the covariant with contravariant indices in the right 
hand side 



lj,\t 



)} 



and the equation is 
where 



(B.160) 



J^ = -i{[4>\D^(f\ - {D 



as in the definition Eq. fH^ . Then the equation represents the component 1 
of the equation of motion, Together with Eq. (lB.119p we have 



where 



i.e. Eq. 



12.5 The Euler-Lagrange equation for the matter fields 

This equation is obtained by functional variation of the action at the matter 
fields, and respectively 0^ 



d 5C 



and 



dx'^ 5 (^) 50 
d 6C 6C 



The matter Lagrangean consists of the kinematical part 

£2 = -tr [{Dy)^ {D,<P) 
and the potential of self-interaction for the matter field 



(B.161) 



(B.162) 



V (0, 0t) 



4^2 



tr 



([[0,0t],0]-t;^0)^([[0,0t],0]-^^ 
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(B.163) 



C,naUer = C2-V{<i>,<P^) (B.164) 

Let us consider the second Euler-Lagrange equation and calculate the func- 
tional derivatives at d^ct)^ . We have 

5C SC2 



d 



5(5o0t) 



^(ao0t) 



tr 



(B.165) 



(B.166) 



Analog calculations give 



S 



and 



SC. 



5(a20t) 



5(5i0t) 
5 



tr 



(L>i0)^ (B.167) 



tr 



(B.168) 



The other term in the Euler-Lagrange calculation imphes the derivatives 

SC 5jC2 SV 



The first term is 

5/12 S 



tr 



£2 = -tr 



(50t 50t (50t 



(B.169) 



(B.170) 



(B.171) 



-tr 

+ 



dx 
J \dy 



+ Ai(j) - (j)A] 
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and we will calculate it in detail. 
S 



_5_ 

_S_ 
1^ 



tr 
tr 
tr 



(B.172) 



'dt 



+ 



at 



ot 



at 



Term by term 



50t 
_S_ 

_S_ 
tr 



tr 



ot: 



dt 



tr [0^0^00] = - (A°Uo0)^ 
tr [-0tA°t0Ao] = {A'^MoY 



5 
5 



dt 



T [ d(j) 
dt 



Summing up the terms 



tr = (^°T ( W 

tr [A'^^A,] = - {A'^f {<pA,f 



5 



tr 



(90 

at 



- {A^^ A,<f)' + {A^^A,y 
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We apply the transpose and factorize 

5 



50t 



tr 



- (L'oc/.)^ (L'o^)] 



T 



or 



50t 



tr 



By analogue calculations we obtain 

S 



and 



Then 



50t 



tr 



tr 



(L>20)^(L>20) 



5£2 



50t 5(/)t 



tr 



- (^00)^ (L'o^) + (A0)^ (A0) 



Further 



5£^_SC2 SV 
~ ~ 
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and the Euler-Lagrange equation results by combining Eqs. flB.166p . flB.167p 
flB.1681) and flRTTl 

d 5C _ 5£ 

Y dxt^ J 



+ 



T 



+ 



D2(P,Al 



6V 



We note that 



d 
dx^ 



dIDo^ 



+ 



The other terms have similar form and we obtain 



13 Appendix C : Derivation of the second 
self-duality equation 



The gauge field equation in terms of ± variables (Dunne 
Let us calculate 

F+_ = d+A_ - + [A+, A_] 
using the space variables (1,2) = {x,y). 

d . d 



dx^ ~'~ dx^ ) ^ ^ "^^ 
+ [Ai + iA2, Ai - iA2] 



(C.l) 



(C.2) 
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dAi Ml 



.dA2 dA2 



dx^ ' dx"^ dx^ dx"^ 
dAi .dAi .dA2 dA2 



— i- 



dx^ dx"^ dx^ 
-I [Ai,A2]+i\A2,Ai] 



dx^ 



-2i 



dA2 dAi 
dx^ dx'^ 



12 



12 



On the other hand, we have the equation of motion Eq.(H6 



2k£°^2Fi2 = if 







from which we derive 



12 



where we can use the second of the Eqs.f lU^ . vahd at self-duahty 



2k 



We calculate separately the terms 



{Do 



2k 



2k 



,2 Ut 
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We can prove by expanding the commutators the equahty of the first terms 
from the curly brackets of right hand sides of the two equations. From 
Eq.(la8]) 



— /At U /Atl 



(C.IO) 



From Eq. flOQl) 



smce 



then 



[0,0t] _ [0^0t] 0t^0] 
0t [0, 



(c.ii) 



We note that the first and the last terms in fIC.lOp and fIC.lip are the same. 
The other terms in fIC.lOp are 



^ 00^ - 0^ 



(C.12) 



and from (IC.lip 



(C.13) 
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and the two expressions are identical. This means that 



' ' 7 ' 



- 7,2 Ut 



and 
and 



(C.14) 
(C.15) 

(C.16) 



where we replace the expression of Dq 



2i 

K 



(C.17) 



K 



14 Appendix D : Notes on definitions 

The information about the algebraic structure invoked in the present model 
can be found in [32] • A simple group is a group that does not have invariant 
subgroups, except of the identity and the whole group. 

A simple algebra is an algebra that does not have proper ideals. 

A semi-simple algebra is an algebra that can be written as a direct 
sum of simple algebras. 

U (1) is not simple. 

The dimension of a simple Lie algebra is the total number of linearly 
independent generators. 

The rank of the algebra, r, is the maximum number of simultaneously 
diagonalisable generators of a simple Lie algebra. 

In the Cartan-Weil analysis the generators are written in a basis where 
they can be devided into two sets: 

• the Cartan subalgebra, which is the maximal Abelian subalgebra of 
G. It contains r diagonalisable generators Hi, i = 1, ...,r 



[H„H,]=0, z,j = l. 



(D.l) 
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• the remaining generators of the algebra G are defined such as they 
satisfy the eigenvalue problems 

[H,,E^]=aiE^, i = l,...,r (D.2) 

It results that the constants can be considered structure constants of 
the algebra in the Cartan-Weil basis. 

For each generator there are r constant numbers, a^, ^ = 1, ...,r ; if 
we consider a space with dimension r, then the set of points (ai,a2, ....ar) 
corresponding to one generator E^ is a point in this space. This space is 
called root space and the name root comes from the fact the the vector 
(ai, 0:2, cxr) is obtained by solving the equation (ID.2p . an eigenvalue prob- 
lem. 

Two problems are connected and are treated together using the Dynkin 
diagrams of the simple algebras: 

1. to classify all possible systems of roots for the algebras of a given 
rank r; 

2. to find all possible irreducible representations of a simple group 
G. This means to identify a system of physical states on which the 
generators E^ are acting (the states belong to a Hilbert space) with 
the property that these states are transformed between them (or, the 
system of states is closed under the action of the generators E^). These 
states are taken as the basis for an irreducible representation. 

Considering the physical states which are the basis of the irreducible rep- 
resentation, I A), they can be labelled by the r eigenvalues of the diagonalisable 
generators 

H,\\) = \\\) , z = l,...,r 
The set A is called the weight of the representation vector. 



15 Appendix E : Expanded form of the first 
equation of motion 

The first equation of motion is 

D,D-4. = 1^ (E.1) 
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As explained by Dunne [oo] the derivative of the potential V is obtained from 
the functional variation to (j)^ and for this we need the expanded form of V. 
The potential is given initially in terms of the trace of the operators 



1 



V{<p,<p^)=^tr[{[[<p,<p^],<P]~v'<py{[[<P,<P^],<P]-v' 
In the equations of motion we will treat separately each term: 

([[0,0t],0]_^20)t([[0,0t],0]_^2^) 



The first term, of sixth degree 

- (f)\ 



t _|_ ^^t^^''' 



t + /A/At 



X 



+< 
+ 



+ 



(E.2) 
(E.3) 



(E.4) 



(E.5) 
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We remark that in Eq. flE.5p the terms five and seven 



tr (-20V0V00"^ + 

terms six and eight 



2tr ^t^t^^^^t^ 



terms second and nine 



third 



four 



tr (-20V0W</> + (E.6) 
_tr f'/At^^t^t^ 



(E.7) 
(E.8) 



can be grouped. Collecting the terms we have 



[[0,0t]^0]^[[0,0t],0] (E.9) 
_tr /'/At^^t^^^t^ 

+2tr ^ 

—"tj- f' /At/A/At/At/ 



-2tr 

-2tr ^'^^t^^^'f' 



The third term, after two permutations of the first two factors, is identical 
to the fifth term 

-tr - 2tr (E.IO) 

= -tr - 2tr f ^t^t^^^t. 

= -3tT ((f)UU(J)(l)h 



The first and the fourth factors are equal 



-tr - 2tr (E.ll) 

— 3tr ^"^^^t^^^'^'i 
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Then from the first, sixth degree product, we obtain 



[[0,0t],0]t[[0^0t]^0] (E.12) 
-3tr f^VV^^V) - 3tr (S^M^Md)^) + 2tr (S^S^MM^) 



The next two terms in the potential (proportional with (— f^)) are ex- 
panded 



— f 



— f ^ 

The last term is unchanged 



2. the linearity for sum of arguments 

tr (A + B) = tr (A) + tr (B) 



(E.13) 



(E.14) 



Now we invoke two properties of the Trace operator: 
1. the symmetry to cyclic permutation 

tr (ABCD) = tr (DABC) = tr (CDAB) = tr {BCD A) (E.15) 



(E.16) 



Then we remark in Eq. flE.13p that the last three terms can be grouped, 
so that the final form for it is 



i^—v'^) tr {40''" ~ '^^''^"l'^^ — ^t^^^t 
i^—v'^^ 4tr Z^"''^^"''' 



(E.17) 
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Adding the contributions to the potential 

4/tV(0,0"f) (E.18) 
= -3tr - 3tr (0^0^00^) + 2tr (0^^000^) 

+ {-v^) 4tr {0V0V - 0W0} 



+vHt 



Consider now the variation of V (0, 0^) to the function 0''' 

This will be calculated by adding a small functional variation to 0^ and 
retaining the first oder: 

perturbed sixth order part (E.20) 

= tr [-3 (0^ + 50^) 0V00V 
-30t (0t + 50t) 000t0 

-30W0 (0^ + 50^) 
-3 (0^ + 50^) 00V00^ 
-30^0 (0^^ + 50^) 000^ 

-30V0V0 (0^ + (^00 + 
+2 (0^ + 50^) 0^000^ 

+20t (0t + 50t) 0000t 
+20t0t000 (0t + 50t)] 

This gives, after permuting the small 50^ to the left 

perturbed sixth degree part (E-21) 
= sixth degree part + 
+tr [-30V00V (^0^) 
-3000V0^ (50^) 
-300VV0 (50^) 

-300'''000''' (5(p^) 
—3000^0^0 ((^0^) 
-30V0V0 {0) 
+ 20^000^^ (50^) 

+20000V^ (50^) 

+ 20W00 (^0^)] 
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This is symbolically written 



perturbed sixth degree part (E.22) 
= sixth degree part 
+tr [A (50t)] 

where A is given in (lE.34p . 
The fourth order part is 

perturbed fourth degree part (E.23) 
= 4tr [(0^ + 50^) 

-0t (0t + 50t) 00] 

or 

perturbed fourth degree part (E.24) 
= fourth degree part + 
+ [-v"^) 4tr [00V (50"^) 

-000t (50t)] 

This can be written 

perturbed fourth degree part (E.25) 
= fourth degree part + 
4tr [B (50^)] 

where 

B = 0(t) + 00V - 0V0 - 000^ (E.26) 

The last part 

perturbed second degree part (E.27) 
= [(0^ + (50^) 0] 
= second degree part 

+t;^tr [0 (50"^)] 
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The three terms have the sum 



V (0, 0+ + (50+) 



+tr [A 

+ (-^^) 4tr [5 
WtT [0 (50+)] 



We introduce a short notation 



C = ^ + 



(- V) S + -uV 



and we have 



V (0, 0t + 50+) = F (0, 0+) + tr [C (50+)] 



V ^21 C^22 '^^k 5<Pl2 ) 
^ Cn50+i + Ci250L Cn50+2 + Ci250: 

^ V ^^215011 + C2250L C^2l50l2 + C^2250i 
= Cn50+i + Cl250^i + C2i50+2 + C2250^2 



A = -30+000+0 - 3000+00+ - 300+0+00 
-300+000+ - 3000+0+0 - 30+00+00 
+20+0000+ + 20000+0+ + 20+0+000 



The last term is 



Prom here we can derive 




In detailed form 



[A+{-Av^)B + v^(f\' 



where 
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and 

B = 2#V - 0t00 - (E.35) 

and the last term is 

(E.36) 
We must make a rearrangement of terms in 

C = ^ + (-4^;2) B + v^cj) 

in order to obtain the last form of the equation of motion. 
NOTE 

For comparison and for an easier analysis we can use the Abelian version 
as a suggestion. The Abelian version of this arrangement is 

(101^-^^) (310^-^^)0 

This may work for example for the second part (proportional with (— f^)) 

(-4^;^) B 

and this is similar with the Abelian version 

Obviously the last term is the same 

v'^cf) <-> {—v^Y (f) 
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